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ABSTRACT 

A digital computer program for analyzing the electromagnetic design of horaopolar 
inductor alternators is presented. The program, which is written in FORTRAN IV pro- 
gramming language, is described in general terms. The calculational methods are either 
outlined briefly or appropriate references are cited. Instructions for using the program 
are given and typical program input and output for a 15-kVA alternator are shown. Cal- 
culated results for this and two (nearly identical) 80-kVA alternators are compared with 
experimental data. In general, considering the many assumptions and approximations 
which are made in the calculational methods, it is felt that reasonable agreement has 
been obtained between the test data and calculated results. 
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DESCRIPTION AND EVALUATION OF DIGITAL-COMPUTER DESIGN-ANALYSIS 
PROGRAM FOR HOMOPOLAR INDUCTOR ALTERNATORS 
by David S. Repas and Gary Bollenbacher 
Lewis Research Center 

SUMMARY 

A digital computer program for analyzing the electromagnetic design of homopolar 
inductor alternators is presented. The program, which is written in FORTRAN IV pro- 
gramming language, is described in general terms. 

The method of calculation is either outlined briefly or appropriate references are 
cited. The items that are calculated by the program include the open- circuit saturation 
curve, the field-current requirement at various loads, losses, efficiency, and react- 
ances. Instructions for using the program are given, and typical program input and out- 
put for a 15-kilovolt-ampere alternator are shown. Calculated results for this and two 
(nearly identical) 80-kilovolt-ampere alternators are compared with experimental data. 
The comparison shows that the maximum difference between calculated and experimental 
data is 7 percent for field currents and 0. 7 percent for efficiency at rated load. 

An alphabetical list of major FORTRAN symbols, the complete program listing in- 
cluding flow charts, and a list of input variables with definitions are given in the appen- 
dixes . 


INTRODUCTION 

The application of the digital computer to the design of alternators has found wide 
acceptance within the electric machinery industry. However, specific computer pro- 
grams that have been written remain for the most part proprietary. 

In 1964 work sponsored by the NASA resulted in a report (ref. 1) that contained 
eight design manuals and eight digital computer programs for analysis of most major 
types of alternators. The programs are written in the FORTRAN II programming lan- 
guage for use on an IBM- 1620 computer equipped with an on-line card reader and a type- 
writer console for input and output. 

These programs suffer from two shortcomings. The first is the limitations imposed 



by the equipment for which it was written. The second and more serious shortcoming is 
that, for most of the programs, accuracy had never been thoroughly verified by compar- 
ing calculated results with experimental data. Both shortcomii^s were remedied for one 
of the eight computer programs. The homopolar inductor program was chosen because 
of the interest in this alternator for use in space-power systems and because of the ready 
availability of experimental data for three different homopolar inductor alternators. 

Elimination of the shortcomings required numerous program modifications. These 
modifications included converting the program to the FORTRAN IV programming language 
for use on an IBM- 7094 computer and rewriting the input and output statements to utilize 
high-speed peripheral equipment. The required input data to the program were substanti- 
ally reduced, and checks for obvious errors in the input data were added. The output was 
clarified to the point of being self-explanatory. 

More significant were the modifications found necessary when results of computer 
calculations were compared with experimental data for the 15-kilovolt-ampere Brayton 
cycle alternator (refs. 2 and 3) and for the two 80-kilovolt-ampere SNAP- 8 alternators 
(refs. 4 and 5). All three of these alternators are rated at 120/208 volts, 400 hertz, and 
12 000 rpm. To obtain satisfactory agreement between experimental and calculated re- 
sults, modifications were made in the magnetic, reactance, and efficiency calculations. 

As shown in this report, the final version of the homopolar inductor alternator com- 
puter program gives calculated results that agree favorably with experimental data for all 
three alternators. The program may be used both for analyzing the electrical design of 
specific alternators and for parametric studies of alternators for auxiliary power gener- 
ating systems. 


CO/VIPUTER PROGRAM DESCRIPTION 
General Description 

The homopolar inductor alternator computer program is an analysis program. This 
means that the program accepts as input a complete electromagnetic alternator design; 
from this, it calculates losses and efficiency, the open-circuit saturation curve, field- 
current requirement at various loads, several reactances, and weights of electromagnetic 
components. The results of the calculations, together with the input, are then printed out 
to provide a complete, self-explanatory record. 

The program may be used with any computer system that accepts FORTRAN IV. 

For program execution, approximately 13 000 storage locations are needed. At the Lewis 
Research Center, the program has been used on the 7044-7094 Mod II direct couple sys- 
tem using a FORTRAN IV, version 13 compiler. For this system, typical pre -execution 
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and execution times for the program are 1. 0 and 0, 04 minute, respectively. 

The computer program consists of a main program and three subroutines. The sub- 
routines were necessary because one long program would have been too large to compile 
with the available core storage locations. 


Description of Alternator to Which Program is Applicable 

The basic alternator configuration for which the computer program was written, 
with each major electromagnetic component identified, is illustrated in figure 1. As 
shown, the alternator consists of two laminated stators separated by a toroidal field coil. 
Surrounding both stators and the field coil is the yoke. The armature winding passes 
through both stators and under the field winding. 

The rotor is constructed with saliences or poles on each end; all north poles are at 
one end and all south poles at the other. As in a conventional salient-pole alternator, 
the centerlines for the north and south poles are 180 electrical degrees apart. 

A number of assumptions, in addition to those implicit in the geometric configura- 


I- Armature winding 



/ / i-Stator 

/ I 

Field coil “* *-Stator back-iron 

or stator core 

Figure 1. - Cutaway view of homopolar Inductor alternator. 
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tton, are made regarding the alternator. These assumptions are 

(1) Shaft, poles, and pole head are made of same magnetic material 

(2) Alternator armature winding is three-phase Y-connected 

(3) Both stators are made of the same material 

(4) Distance between stators is the same as the field coil width 

(5) Field coil is confined to the toroidal space bordered by a stator on each side, the 
yoke on the outside, and the armature winding on the inside 

(6) Alternator has only one field winding. 

In contrast to the restrictions imposed on the alternator by the preceding assumptions, 
there are several options that are available to the program user. These options, which 
increase the applicability of the program, are 

(1) Armature conductors may be round or rectangular 

(2) Field conductors may be round or rectangular 

(3) Armature conductors may consist of any number of strands 

(4) Yoke, rotor, and stator may each be made of a different magnetic material 

(5) Damper windings may or may not be present 

(6) If damper windings are present, the damper bars may be either round or rec- 
tangular 

(7) Five different slot configurations may be used 

(8) Three different yoke geometries may be used. 


/Vlethod of Calculation 

This section of the report will outline in general terms the method of calculation used 
in the computer program. However, due to the length of the program and the large num- 
ber of equations involved, specific equations will not, except in a few instances, be given. 
Instead, references for the major design analysis equations are given. Reference 1 is 
particularly applicable. 

More detailed information and specific equations may be found in the program list- 
ing in appendix A . To assist in locating specific information in the listing, COMMENT 
cards are used freely to identify the major calculations. Of further value is appendix B, 
which is an alphabetcal listing of the major FORTRAN variables including definitions and 
units, and the flow charts for the main program and two of the three subroutines included 
in appendix A. 

Magnetic calculations . - A cross-sectional view of a homopolar inductor alternator 
is given in figure 2. For clarity, a two-pole alternator is shown. The main flux path in 
the alternator is shown by the solid arrows, and the leakage flux paths are indicated by 
the broken arrows. An additional leakage flux from the rotor to the stator between 
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Figure 2, - Homopolar-inductor-alternator basic configuration 
and flux paths. Leakage flux across field coil, (p-^; leakage 
flux from stator to stator, (p 2 \ leakage flux from stator to rotor 
end extension, (p^^. 


the rotor poles is also present. This path is shown in figure 3, which is a developed end 
view of the alternator. 

The main flux flows from a rotor north pole, across the air gap and then radially 
through the stator teeth and stator back iron. It then goes axially throi^h the yoke to the 
other stator stack where the flux path is completed through the stator laminations, air 
gap, and rotor. 


Stator 



Figure 3. - Homopolar-inductor-alternator end view 
showing leakage flux between poles. Leakage flux 
between poles from rotor to stator, (p^. 
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Field coil 



Figured. - Equivalent magnetic circuit for homopolar inductor 
alternator at no-load. 


An equivalent magnetic circuit for the homopolar inductor alternator is given in fig- 
ure 4. The various leakage fluxes and permeances which are considered in this program 
are shown. In this report, the laminated stator back iron is referred to as the stator 
core. The rotor shaft is the cylindrical part of the rotor and excludes the poles (fig. 1). 

Some of the more important equations and assumptions used to determine field cur- 
rents for various load conditions will be described in this section of the report. The 
complete equations for the magnetics calculations can be found in the FORTRAN program 
listing for subroutine MAGNET which is included in appendix A. 

The method of calculation used to determine the field current at no-load will be de- 
scribed first. The flux distribution in the air gap at no-load is shown in figure 5. In the 



Figure 5. - No-load rated-voltage airgap flux distribution for 
homopolar inductor alternator. 
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following discussion, the useful flux in the air gap and poles is taken to be the flux that is 
present, excluding the leakage flux between the poles from the rotor to the stator (p^. 

In determining the useful flux in the air gap and poles, a hypothetical total flxix cp^ 
is first calculated. This hypothetical total flux is assumed to have a constant flux density 
over the entire pole pitch; that is, the shape of the field form is assumed to be rectan- 
gular (ref. 6). 

From the equation for the induced voltage in the armature winding of a synchronous 
machine, is calculated. This calculation takes into account the fact that the winding 
is pitched and distributed and that the actual flux wave is not a true sinusoid. The flux 
density in the air gap B due to the useful flux is 

o 


TTld 


where 

I length of one stator stack 
d inside diameter of the stator laminations 

The air gap magnetomotive-force drop from is then 


F 


g 



^o 


where 

gg effective length of the air gap 
permeability of air 

The useful flux per pole cp is 

Jb^ 


(p 


P 




P 


C 


P 
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where 


P number of poles 

Cp ratio of the average to the maximum value of the field form (ref. 7) 

From (p^, the flux densities and magnetomotive -force drops, due to the useful flux 
in both the poles and stator teeth, are determined. The effect of <p the leakage flux 
between poles from the rotor to the stator, on the air gap, pole, and teeth flux densities 
and magnetomotive -force drops must now be included. 

It is assumed that this leakage flux density is constant around the stator periphery 
(fig. 5). Also, is the product of the sum of the air gap plus pole plus teeth 
magneto motive- force drops and the permeance of the leakage path. The effect of cp^ is 
to increase the flux densities and, thus, the magnetomotive-force drops in the air gap, 
poles, and teeth. Since the magnitude of and these magnetomotive- force drops are 
interrelated, an iteration process is involved in determining cp^. 

Once the preceeding part of the magnetics calculations is completed, the rest of the 
procedure is fairly straightforward. A flow chart that gives the order of the entire mag- 
netics calculations is given in appendix A. The magnetomotive-force drops in the mag- 
netic parts of the alternator are determined in the program from the material magnetiza- 
tion curves. These curves are an input to this program. 

The magnetics calculations are also made for several alternator loads at the load 
power factor specified in the program input. Rated terminal voltage for the alternator is 
assumed for these calculations. For load conditions, some modifications to the no-load 
calculation method must be made. As shown in reference 1, the air-gap magnetomotive 
force under load F^^ will increase from the no-load rated-voltage value. Neglecting 
the effect of 

^ m’ 


F_, = e j • F_ 
gl a g 


where 


e^j = WX^ sin ij/ + cos(ij/ - 0) 

. sin 6 + WX 
rf/ = tan" ^ ^ 
cos 0 


and 

W load at which Fg^ is to be calculated, per unit 
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direct-axis synchronous reactance, per unit 
e cos“^ (power factor) 

Xq quadrature -axis synchronous reactance, per unit 

Also, from reference 1, the flux per pole under load will increase from the no-load, 
rated-voltage value. Again, neglecting cp^, the flux per pole imder load <ppi is 

^ * *^p 


where 


gx = ®d ■ ^^ad ^ 


where X^^ is the direct -axis armature reaction reactance. Now, cp^ is a function of 
the air gap, pole, and teeth magnetomotive -force drops and of the demagnetizing 
magnetomotive -force due to the armature current. 

Using these modifications, the magnetic characteristics of the alternator for load 
conditions can now be determined. The procedure is essentially the same as presented 
for the no-load case. 

Efficiency and loss calculations . - Individual losses and efficiency are calculated at 
several loads of increasing magnitude, continuing until the alternator saturates or vmtil 
calculations have been completed for five loads. While the first load at which loss cal- 
culations are made must always be zero per unit, the program user has the option of 
specifying any or all of the remainii^ four loads. These loads are designated by G 
within the program (G is in per unit). 

Rated voltage and power factor, as defined by the program input data, are assumed 
throughout the loss and efficiency calculations. The individual losses, that are calculated 
by the program, aloi^ with the method of calculation or references, are listed below. 

Field conductor losses and armature conductor losses: These losses are given by 

9 

the expression I R where I is the dc or rms current in the winding, as appropriate, 
and R is the dc winding resistance corrected for the winding temperature. Correcting 
the winding resistance for temperature involves several assumptions: 

(1) The average no-load winding temperature known or can be estimated. 

(2) The average rated-load winding temperature Tj^l known or can be estimated. 

(3) The average winding temperature is a parabolic function of the current in the 
winding. 

With these assumptions, the winding temperature Tq at any load G is 
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where 




+ T 


NL 


Ij^L current at rated load 

Inl current at no-load, equal to zero for armature winding 
Iq current in winding at load G 

For the armature winding is, of course, zero. (If in the program, 1.0 per unit 

load (G = 1.0) is not one of the loads for which losses are calculated, then the above 
equation is only approximately applied to the field temperature calculations. ) 

Eddy losses: References 1 and 8 present discussions of armature conductor eddy 
losses. 

Pole-face losses: For no-load pole-face-loss calculations, see references 1 and 9; 
for pole-face-loss calculations at any other load see reference 10 (eq. 22), 

Damper losses: No-load damper losses are calculated as shown in reference 11 
using the "cold” damper-bar temperature; for damper bar loss calculations under load 
(ref. 10, eq. 22), the "hot" damper bar temperature is used regardless of the magnitude 
of the load. The cold and hot damper-bar temperatures are inputs to the program. 

Stator core loss and stator tooth loss: The respective equations used to calculate 
these losses are 

Stator core loss = k(Stator core weight) density 

\ BK 

Stator tooth loss = k(Stator tooth weight) (WL)/-!^^^^ 


where 

k empirical constant equal to 3. 0. (This constant is variously stated in the literature 
to range from 1. 5 to 3. 0. The 3. 0 value was chosen in this program because it 
provided the closest agreement between experimental and calculated values. ) 

WL core loss at flux density BK and at rated alternator frequency, W/lb 

BK flux density at which WL is measured 
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and where weights are given in pounds. 

Windage loss: If an accurate value for windage loss is known, it may be read into 
the program for use in the efficiency calculation. If the windage loss is not read into the 
program, it will be assumed to be zero. The program user may also elect to have the 
program calculate an approximate value for windage loss. In that case, the equation used 
(ref. 1) is 

W = 2. 52X10" ®(d^- ^n^' h) 


where 

W windage loss, W 
d rotor diameter, in. 
n rotor speed, rpm 
I pole length, in. 

This equation assumes that the gas surrounding the rotor is air at standard pressure and 
temperature. For gases other than air at standard pressure and temperature, windage 
losses may be calculated by the method given in reference 12. 

Miscellaneous load losses: These losses are assumed to be 1 percent of the kilovolt- 
ampere output of the alternator at load point G. 

Efficiency: At each load efficiency is calculated from 

Efficiency = - Alternato r power output ^ mn 

Alternator power output Losses 

where both alternator power output and losses are expressed in watts. 

Reactances . - In the program, the following reactances are calculated; 

(1) Armature winding leakage ^al 

(2) Direct-axis armature reaction Xad 

(3) Quadrature- axis armature reaction 

(4) Direct- axis synchronous 

(5) Quadrature- axis synchronous X^ 

(6) Field leakage Xj 

(7) Direct-axis transient X^ 

The armature winding leakage reactance is the sum of the slot leakage and end wind- 
ing leakage reactances. The slot leakage reactance is determined from formulas given 
in reference 7, but the end winding reactance is calculated using the method of refer- 
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ence 13. Both the direct and quadrature-axis armature reaction reactances are deter- 
mined from the method given in reference 7. The synchronous reactances are determined 
in the usual manner; that is, + ^al and Xq = Xg^q + Xg^j. 

The field leakage reactance is determined from the permeances of the alternator 
leakage paths. These paths are shown in figures 2 to 4. The field leakage permeance 
is 


^ =- 

f 


1 


m 


where 


permeance of leakage path across field coil 
permeance of leakage path from stator to stator 
^2 permeance of leakage path from stator to rotor and extension 
P number of poles 

"^m pe^^J^62.nce of leakage path between poles 
The field leakage inductance Lj is then 


Lj - Nj • 


where 

Nj number of field turns 

The field leakage reactance referred to the field winding Xjj is 

Xjj = 27rf • Lj 


where 


f rated output frequency of alternator 

The field leakage reactance referred to the armature is then 



where 
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where 


N = ® 


N„ • N„ 




2 • M • C 


effective armature winding turns 
N number of slots 

S 

Ng conductors per slot 
kp pitch factor 

distribution factor 
M number of phases 

C number of parallel circuits 

The direct-axis transient reactance is calculated by the usual method 




• X. 


ad 




Skew factor calculation . - The usual skew factor formula for conventional alternators 
having only one stator stack does not apply to a homopolar inductor alternator. A new 
equation, which takes into account the stator stack separation, had to be derived for use 
in this computer program: 


Skew factor = 



ST! 

■ 

sin ° 


2T„ 


P 





where 

Tp . 

Sq stator slot skew measured at the stator bore (for one stator stack) 
b distance between two stator stacks 
length of one stator stack 

The preceding equation reduces to the usual formula when the stator separation is zero 

(b = 0) providing that it is recognized that setting b = 0 gives a stator stack of length 

2Z and a total slot skew of 2s . 
o o 
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HOW TO USE COMPUTER PROGRAM 
Input Data Requirements 

To use this computer program for the analysis of a homopolar inductor alternator the 
complete electromagnetic design of the alternator must be known. This includes physical 
dimensions, armature and field winding parameters and the magnetic characteristics of 
the materials to be used in the stator, rotor, and yoke. The design information must 
then be transferred onto data cards for use with the program. A typical set of data cards 
is shown in figure 6. It consists of three material decks. The material decks must be in 
the order shown in the figure, that is, stator material, rotor material, and yoke mate- 
rial. There must be exactly three material decks in each data deck even if two or all 
three materials are identical. 

If more than one alternator design deck is included in the data deck, the program 
will treat each design deck independently. Each will result in a separate alternator analy- 
sis complete with an individual output record. However, the same material decks will be 
assumed to apply to each alternator design deck. 


Any number of 
alternator 
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Preparation of Material Decks 


A material deck consists of five cards. The first card contains the material name. 
This serves two functions: it identifies the material deck, and it is read by the computer 
and stored for later printout on the output record. The remaining four cards contain in- 
formation about the magnetization curve of the material specified on the first card. This 
information allows the approximate reconstruction of the magnetization curve durii^ pro- 
gram execution. Table I summarizes the information pertaining to each data card of a 
material deck. 


TABLE I, - FORMAT AND TYPE OF DATA REQUIRED ON 
MATERIAL DECK DATA CARDS 


Card 

Format 

Information contained on card 

1 

6A6 

Material name 

2-5 

8F10. 1 

Coordinates from material magnetiza- 



tion curve 


To illustrate preparation of a material deck, AISI 4620 steel (hardened) will be used 
as an example. The first card of this material deck will appear as shown in figure 7, 

The material name should start in column 1 and may extend up to colunm 36, 

To prepare the remaining four cards of the material deck, the magnetization curve 
of the material is needed. The magnetization curve for AISI 4620 steel (hardened) is 
shown in figure 8. The units must be kilolines per square inch for the magnetic flux 
density and ampere-turns per inch for the magnetizii^ force. Fourteen points on the 
curve must then be chosen. In the figure, 13 points are indicated by data symbols; the 
14th point is the origin. These points are listed in the table insert. Careful attention must 
be paid to the sequence in which the numbers are punched onto data cards. The first num- 
ber must be the maximum flux density of the points chosen. In the example, this value is 
128 kilolines per square inch. This is followed in ascending order, by alternate values 
of magnetic flux density and magnetizing force. Again, in the example, with reference 
to the table Insert, the values appear in the following sequence on the data cards: 128, 

0, 0, 2, 5, 5, 10, . . . 110, 115, 128, 300. llie complete material deck for AISI 4620 
steel (hardened) is shown in figure 7. 

During program execution, the original magnetization curve is approximately recon- 
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Figure 7. - Material deck for AISI 4620 steel (hardened). 
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structed by interpolation between points. The interpolation assumes a straight line on 
semi-log paper between data points. 


Preparation of Alternator Design Deck 

The alternator design deck contains ail the dimensions, the geometric configuration 
(in numerical code), and the winding parameters needed for an electromagnetic analysis 
of the alternator design. Unlike the material decks, which are read according to a 
FORMAT statement, the alternator design decks are read with a READ statement refer- 
encing a NAMELIST name. For each NAMELIST name one or more data cards are re- 
quired to numerically define the variables included in that NAMELIST name. In all there 
are 11 NAMELIST names. Each name is suggestive of the type of variables Included in 
its list. Table II lists the NAMELIST names in the order in which they must appear in 
the alternator design deck and indicates the t5?pe of information conveyed by the variables 
belonging to that NAMELIST name. Detailed information about each NAMELIST name is 
provided in appendix C . 

Preparation of an alternator design deck will now be illustrated with the construction 
of a t 3 rpical data card for the NAMELIST name DAMPER. Ibe data that will be used is 


TABLE n. - SUMMARY OF NAMELIST NAMES USED IN 


ALTERNATOR DESIGN DECK^ 


NAMEUST 

Type of information included 

b 

name 


RATING 

Rated kVA, power factor, voltage, rpm, etc. 

STATOR 

All stator dimensions but not including slot 


dimensions 

SliOTS 

Specifies type of slot and slot dimensions 

WINDNG 

Fully describes armature winding 

AIRGAP 

Gives air gap dimensions 

CONST 

Gives various constants needed for internal 


calculations 

ROTOR 

Gives pole and pole head dimensions but not 


including damper winding 

DAMPER 

All variables concerning damper windings 

SHAFT 

All shaft dimensions 

YOKE 

Yoke dimensions and type of yoke 

FIELD 

Includes all field coil parameters 


^For detailed information, see appendix C (table vn). 
^Presented in the order in which they must appear in die 


alternator design deck. 
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0.070 






Damper bar length 

2.12 

Damper bar pitch 

0.32 

Number of damper bars per pole 

9 

Damper bar material 

Copper 


(a) Damper bar and damper bar slot design for Brayton -cycle alternator. 

(All dimensions are in Inches.) 

BN 9 

WO 0.030 

HD 0.070 

DD 0. 100 

H Not read in when damper bars are round 

B Not read in when damper bars are round 

SB 2.120 

TB 0.320 

T33 20° C is acceptable since true temperature is (by assumption) unknown 

T3 130 

RE Not read in since 0. 694 is sufficiently accurate 
ALPHAE Not read in since 3. 93x10'^ is sufficiently accurate 



WD=0.03» HD=0.07» DD=0.10f BN=9» SB=2. 12, TB=.32» T3=130 $ 


FORTRAN STATEMENT 

jam II I 0 Qooi I I aoia oi a'oaia a eooaoa a oeoaa o ooaa o o o o o ooo o o ajiA^ 
h 11 1 1 1 n 1 1 1 1 1 1 1 1 1 n 1 1 1 1 11 n i n 1 1 1 1 1 1 1 1 1 1 1 11111111111111 inii** 

' ■» j j 2 2 2 2 : 2 2 2 2 2 ? J > ♦ • ' • 


(b) Numerical values of DAMPER variables and appearance of data card. (See 
table Vll(h) for definitions of FORTRAN symbols. ) 

Figure 9. - Preparation of data card for NAMELIST name DAMPER. 


for the 400-hertz, 15 -kilovolt -ampere, 120/208-volt Brayton cycle alternator (refs. 2 
and 3). Figure 9(a) gives all pertinent design data for the Brayton cycle alternator dam- 
per circuit. Figure 9(b) shows how the design data are related to the variables of 
NAMELIST name DAMPER (table Vn(h), appendix C) and how these data are transferred 
to the data card DAMPER. 

Data cards for the remaining NAMELIST names are prepared in a similar manner. 
To illustrate the result, a complete data deck listing for the Brayton cycle alternator fol- 
lows. 
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II 


SDATA 

SILICON STEEL (.007 IN. LAMINATION) 

129. 0. 0. 3.3 .4 12.9 .8 23.9 

1.2 36.1 1.6 45.1 2. 61.3 3. 68.4 

4. 77.5 6.1 80. 8.1 93.5 60.6 103.2 

181.8 109.8 303. 129. 707. 

AISI 4620 STEEL (HARDENED) 

128. 0. 0. 2. 5. 5. lO. lO. 

15. 20. 21. 30. 24. 50. 28. 60* 

32. 70. 36. 80. 43. 90. 55. 100. 

75. 110. 115. 128. 300. 

INGOT IRON 

125. 0. 0. 2. 1.5 6. 1.8 8.5 

2. 30.5 2.6 48. 3.3 62.5 4,2 75.5 

5.7 59. 9.2 97. 14. 104.5 26. 114. 

98, 121. 210. 125. 300. 

SRATING VA=15» EE = 208, F = 400» IPX=4, PF = 0.8, G=0 , . 5 . 1 . . 1 . 5 , 2 . $ 

SSTATOR DI=5.28, DU=8.68. CL=2.00. LTS=.007. WL=8.6. BK=77.4. SF=0.90 S 
SSLOTS ZZ=2, B0=.065» BS=.171, H0=.04, HX=.482, HS=.62, HT=.035» HY=0.127» 
100=48 $ 

SWINDNG RF=1. SC=8. YY=8» C=2. DW=.140. SN=2. SN1=2. DW1=.0250» CE=.12» 
SD=.0290. Tl=114.5. Tll=93.5 S 
SAIRGAP GC=.040 $ 

SCONST $ 

SROTOR PL=1.88, HP=0.85, HP1=1.0» PE=.700, BP=2.37. LTR1=0.014 $ 

SDAMPER W0=0.03. HD=0.07. DD=0.10. BN=9» SB=2.12. TB=.32, T3=130 S 
SSHAFT DSH=3.53» DISH1=1.00. ALH=2.28 $ 

SYOKE TYPY=1. TY=.44 $ 

SFIELD PC0IL=6.56» DC0IL=8.18» PT=515. RD=.0571, T2=113, T22=l00.5 S 


Typical Computer Program Output 


In this section, the output, which resulted from the input data shown in the proceeding 
section, is presented. This output is typical, although the actual program output format 
will vary somewhat, depending, for example, on the type of slot or yoke configuration 
specified in the input data. 


**H0M0PCLAR INDUCTOR ALTERNATOR** 


ALTERNATOR RATING 


ALTERNATOR KVA 15.0 

LINE-LINE VOLTAGE 208. 

LINE-NEUT. VOLTAGE 120. 

PHASE CURRENT 41.64 

POhER FACTOR 0.80 

PHASES 3 

FRECUENCY 400. 

POLES 4 

RPM 12CG0.0 
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STATCR SLOTS 


TVPE-PARTIALLY CLOSED 


-BO- 





HO 


« 

» 




— — 

— 

— * 

» 

BO 

0.065 

INCHES 



• 

• 

BS 

0.171 


HT 


* 

* 

HO 

0.040 



« 


• 

HX 

0.482 


— 

-* 



HT 

0.035 



• 



HW 

0.052 


HW 

« 



HS 

0.620 


— 

-» 

*••• 

*««• 





» 

» 

• 

NO. OF SLOTS 

48 



• 

• 

• 





• 

*»«• 

«••• 

SLOT PITCH 

0.346 

INCHES 

HX 

• 







• 

»•«««•#« « 

SLOT PITCH 




« 

* 

« 4 

AT 1/3 DIST. 

0.373 

INCHES 


• 

« 

• 4 




1 1 

X BS 1 

1 1 


AIR GAP 

MIMMtP AIR GAP C.C40 INCHES 

HAXIRiJP AIR GAP 0.C40 

EFFECTIVE AIR GAP 0.043 

CARTER COEFFICIENT 

STATOR 1.057 

ROTOR 1.014 


ARPATliRE WINDING I Y-CCNNECTED , FORM WOUND) 

STRAND DIMENSIONS 0.1400 X 0.0250 INCHES 

UNINSULATED STRANC HEIGHT (RACIAL) 0.0250 

DISTANCE BTWN CL OF STRANDS (RADIAL) 0.0290 


STRANDS/CONDUCTOR IN RADIAL DIR. 
TOTAL STRANCS/CONCLCTOR 
CONDUCTOR AREA 

CURRENT DENSITY AT FULL LOAD 

COIL EXTENSION BEYCNC CORE 
MEAN LENGTH OF 1/2 TURN 
END TURN LENGTH 
STATOR SLOT SKEW (PER STATOR) 

RESISTIVITY AT 20 DEG. C 

STATOR RESISTANCE AT 25. DEG. C 


2 . 

2 . 

0.0070 SQ-IN. 

2973.99 AMP/SQ-IN. 

0.120 INCHES 
12.030 
5.750 

0 . 

0.6940 MICRO OHM INCHES 
0.0389 OHMS 


NO. OF EFFECTIVE SERIES TURNS 

SLOTS SPANNED 

SLOTS PER POLE PER PFASE 

CONDUCTCRS/SLCT 

NO. OF PARALLEL CIRCUITS 

PHASE BELT ANGLE 


26.54 

8 . 

4.00 

8 . 

2 . 

60 * 


DEGREES 


SKEW FACTOR 
DISTRIBUTION FACTCR 
PITCH FACTCR 


l.OCO 

0.958 

0.866 
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FIELD ^1^DI^G 


CONDUCTOR CI/i^'ETER 
CONDUCTOR AREA 


0.057X INCHES 
0.0026 SQ-IN. 


NO. OF TURNS 

MEAN LENGTH OF TURN 


515. 

23.154 INCHES 


RESISTIVITY AT 20 DEG. C 
FIELD RESISTANCE AT 25. DEG. C 


0.6940 MICRO OHM INCHES 
3.2951 OHMS 


COIL INSIDE DIAMETER 
COIL OUTSIDE DIAMETER 
COIL V^IDTH 


6.560 INCHES 
8.160 
2.280 


STATCR 


STATOR INSIDE DIAMETER 
STATOR OUTSIDE DIAMETER 
OVERALL CORE LENGTH (ONE 
EFFECTIVE CORE LENGTH 
DEPTH BELCW SLOT 

STACKING FACTOR 



5.28 INCHES 


8.68 

STACK) 

2.00 


I.6C 


1.08 


0.90 


NO. OF CCCLING DUCTS 
WIDTH CF DUCTS 


C. 

0. INCHES 


CORE LOSS AT 11 .^ K ILCLINES/SQ. IN. 
LAMINATICN THICKNESS 


8.6 WATTS/LB. 
0.C07 IN. 


ROTOR 


POLE 

eCDY WIDTH 

2.370 

INCHES 


AXIAL LENGTH 

1.880 



STACKING FACTOR 

1.000 


POLE 

HEAD WIDTH 

2.717 

INCHES 


AXIAL LENGTH 

1.880 



STACKING FACTOR 

0.911 



LAMINATICN THICKNESS 

0.014 

INCHES 

POLE 

EMBRACE 

0.700 


POLE 

HEIGHT 

0.850 

INCHES 

POLE 

HEIGHT (EFF.) 

1.000 


ROTOR 

DIAMETER 

5.200 


PERIPHERAL SPEED 

16349. 

FEET/MIN. 

SPEC. 

TANGENTIAL FORCE 

1.076 

LBS/SQ.IN 
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SHAFT 


DIAMETER (ONCER FIELD COIL) 

3.530 INCHES 

INSIDE DIAMETER (CF FOLLCW 

SHAFT) 0. 

DIAMETER (UNDER END TURNS) 

1.000 

LENGTH (BTW. PCLES) 

2.280 

DAMPER BARS (ROUND) 

DAMPER BAR DIAMETER 

0.100 INCHES 

SLCT OPENING kIDTF 

0.030 

SLOT OPENING HEIGHT 

0.070 

DAMPER BAR LENGTH 

2.120 

DAMPER BAR PITCH 

0.320 

NO. OF DAMPER BARS/PCLE 

9 

RESISTIVITY AT 20 DEG. C 

0.694 MICRO-OHM INCHES 


YCKE (TYPE 1) 


1 

1 

1 

1 


YOKE * TY = 0.44 IN. 

• « 



« • • 1 

STATOR • FIELD * STATOR • 1 

« COIL • • 1 

« * ♦ 1 


• • « 


INSIDE YCKE CIANETER 

8.680 

INCHES 

STATOR SEPARATION 

2.280 

INCHES 

kEIGhTS 

STATOR CCND. 

10.380 

POUNDS 

FIELD CCND. 

9.802 


STATOR IRON 

33.092 


ROTCR 

21.966 


YOKE 

22.405 


TOTAL 

(ELECTRCPAGNETIC ) 

97.6A5 



22 



CONSTANTS 


Cl, FUNCANENTAL/MAX. OF FIELC FLUX 1.128 

CP, PCLE CONSTANT 0.711 

CM, DEKAGNETIZATICN FACTOR 0.8A4 

CQ, CRCSS MAGNETIZATION FACTOR 0.5C2 

01, POLE FACE LOSS FACTOR 1.170 


PERMEANCES ( LI NES/ AMPERE TLRN) 

AIR GAP 196.370 PER INCH OF CORE LENGTH 

WINDING LEAKAGE ^ STATOR SLOT 5.078 

STATOR END 9.949 

LEAKAGE 

PM, FROM RCTOR TO STATOR 
(BTWN. ROTOR TEETH) 32.571 

P5, ACROSS FIELC COIL 17.543 

P6, FROM STATOR TC STATOR 27.910 

P7, STATOR TO SHAFT END 6.951 


REACTANCES 

AMPERE CCNDUCTORS/INCH 417.343 

REACTANCE FACTOR 0.586 


STATOR WINCING LEAKAGE 
ARM. REACTION (DIRECT) 
ARM. REACTION (QUAD. ) 
SYNCHRONOUS (CIRECT) 
SYNCHRCNCUS (CUAD.) 
FIELD LEAKAGE 
TRANSIENT 


11.782 PERCENT 
109.504 
57.821 
121.286 
69.603 
75.040 
56.309 


FIELD SELF INDUCTANCE 


1.379 HENRIES 


OPEN CIRCUIT TIME CONSTANT 

(FIELD ONLY) 0.41859 SECONDS 

SHORT CIRCUIT AMPERE-TURNS 1394.682 

SHORT CIRCUIT RATIO 1.C43 

STATCR MATERIAL - SILICON STEEL (.007 IN. LAMINATION) 

ROTOR MATERIAL — AISI 462C STEEL (HARDENED) 

YOKE MATERIAL INGOT IRCN 
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MAGNETIZATION CHARACTERISTICS 
<NG LOADt RATED VOLTAGE) 


TOTAL USEFUL FLUX 

1418.98 

KILOLINES 

USEFUL FLUX/PCLE 

252.19 


FLUX DENSITIES 

AIRGAP (INCL. PML) 

42.77 

KL/SQ-IN 

POLE 

59.12 


TOOTH 

86.38 


CORE 

42.56 


SHAFT (UNDER FLC.) 

58.56 


YOKE (OVER FLD.) 

50.65 


AMPERE-TURNS 

AIRGAP 

6C0.99 

PER STATOR 

POLE 

27.02 


TOOTH 

13.06 


CORE 

2.03 


SHAFT (UNDER PCLE ) 

42.14 


SHAFT (UNDER FLC.) 

71.94 


YOKE 

12.49 



TOTAL 


1454.92 


ALTERNATCR LOAD CHARACTERISTICS (RATED VCtTAGE, 0.80 POWER FACTOR) 


PERCENT LOAD 

0. 

50. 

100. 

150. 

200. 

LEAKAGE FLUX (PPL) 

20.88 

41.09 

64.08 

89.31 

115.20 

AIR-GAP AMPERE TURNS 

600.99 

877.98 

1208.45 

1563.37 

1929.39 

FLUX DENSITIES (KL/SQ-IN) 

POLE 

59.12 

64.33 

69.79 

75.99 

82.74 

TEETH 

86.38 

93.99 

101.97 

111.02 

120.89 

SHAFT (UNDER FLD. ) 

58.56 

67.64 

77.59 

88.72 

100.50 

CORE 

42.56 

46.91 

51.66 

57.10 

63.00 

YOKE (OVER FLC.) 

50.65 

59.96 

70.59 

82.68 

95.54 

TOTAL AMPERE TURNS 

1454.92 

2084.C3 

2893.04 

3875.65 

4913-03 

FIELD CURRENT (AMPS) 

2.83 

4.05 

5.62 

7.53 

S.54 

CURRENT DENS. (FIELD) 

1103.24 

1580.28 

2193.73 

2938.83 

3725.45 

FIELD VCLTS 

12.02 

17.34 

24.79 

35.40 

49.34 

TEMPERATURES (DEG.C) 

FIELD 

100.50 

1C2.89 

113.00 

135.92 

172.78 

ARMATURE 

93.50 

98.75 

114.50 

140.75 

117.50 

RESISTANCES (OHMS) 

FIELD 

4.25 

4.28 

4.41 

4.70 

5.17 

ARMATURE 

C.0492 

0.0500 

0.0523 

0.0563 

0.0618 

EDDY FACTOR 

1.02 

1.02 

l.Ol 

1.01 

1.01 

ALTERNATOR LOSSES (HATTS) 

FIELD 

33.95 

70.16 

139.25 

266.39 

470.69 

WINDAGE 

0. 

0. 

0. 

0. 

C. 

STATCR TOOTH 

98.26 

116.32 

136.92 

162.31 

192.43 

STATOR CORE 

102.55 

124.61 

151.12 

184.60 

224.76 

POLE FACE 

84.97 

90.18 

105.79 

131.81 

168.23 

CAMPER 

0.22 

0.28 

0.32 

0.40 

0.52 

STATCR CCPPER 

0. 

64.98 

272.19 

658.51 

1285.34 

EDDY 

0. 

0.34 

1.30 

2.71 

4.39 

MISC. LCAD 

0. 

75.00 

150.00 

225.00 

300.00 

TOTAL 

319.96 

541.86 

956. 9C 

1631.73 

2646.36 

ALTERNATOR OUTPUT (KVA) 

0. 

7.50 

15.00 

22.50 

30.00 

ALTERNATOR OUTPUT (KW) 

0. 

6.00 

12. OC 

18.00 

24. CO 

ALTERNATOR INPUT (KW) 

0.32 

6.54 

12.96 

19.63 

26.65 

PERCENT LOSSES 

100.00 

8.28 

7.39 

8.31 

9.93 

PERCENT EFFICIENCY 

0.00 

91.72 

92.61 

91.69 

90.07 
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NC-LOAD SATURATION DATA 


VOLTAGE 


PERCENT 

8C.0C 

90. CO 

100.00 

110. CC 

120.00 

130.00 

14C.C0 

145.00 

^0. 

a« 

0. 

LINE-NEUTRAL 

96.01 

1C8.08 

120.09 

132.10 

144.11 

156.12 

168.12 

174.13 

C. 

0. 

LINE-LINE 

166. AC 

187.20 

208.00 

228.80 

249.60 

270.40 

291.20 

301.60 

0. 

0. 

FIELD CURRENT 

2.27 

2.53 

2.83 

3.22 

3.81 

4.50 

5.31 

5.79 

0. 

0. 

FLt> CENS.(KL/SC-IN) 

PCLE 

A7.27 

53.18 

59.12 

65.16 

71.39 

77.71 

84.09 

87.38 

0. 

0. 

TCCTH 

69.07 

77.70 

86.38 

95.21 

104.30 

113.54 

122.87 

127.67 

0. 

0. 

SHAFT 

A6.7C 

52.62 

58.56 

6A.76 

71.49 

78.48 

85.65 

89.51 

c. 

0. 

CCRE 

3A.02 

38.27 

A2.56 

A6.97 

51.63 

56.42 

61.31 

63.88 

0. 

0. 

VCKE 

A0.A7 

A5.A9 

50.65 

56.20 

62.50 

69.20 

76.19 

80.06 

0. 

0. 

AMPERE-TLRNS 

AIRGAP 

A80.53 

5A0.58 

600.99 

662.39 

725. 72 

790.03 

055.02 

886.51 

0. 

0. 

PCLE 

23. A2 

2A.96 

27.02 

29.06 

31.53 

35.29 

40.66 

44.09 

0. 

0. 

TOOTH 

2.56 

3.88 

13.06 

45.87 

123.72 

223.49 

337.56 

417.45 

0. 

0. 

CCRE 

1.65 

1.82 

2.03 

2.27 

2.55 

2.87 

3.25 

3.60 

0. 

0. 

SHAFT 

139.98 

1A7.13 

156.22 

166.48 

178.87 

195.56 

219.19 

234.29 

0. 

0. 

YOKE 

10.78 

11.55 

12. A9 

13.70 

15.21 

17.81 

40.35 

42.12 

0. 

0. 

TOTAL 

1167. 1C 

1301.17 

1A5A.92 

1659.37 

1961.14 

2316.75 

2732.51 

2983.74 

0. 

0. 


^All zeros in a column indicate that some section of the alternator has saturated. Examination of the previous column will generally identify 
which part of the magnetic circuit saturated. 



EVALUATION OF COMPUTER PROGRAM 


Agreement between results of the computer calculations and experimental data was 
determined for three homopolar inductor alternators. These three alternators were the 
400-hertz Brayton cycle alternator and both the preprototype and prototype SNAP- 8 ma- 
chines. A more detailed description of these alternators is given in the following section 
of the report. Test data for the Brayton cycle alternator were obtained from reference 3, 
For the SNAP- 8 alternators, test data were taken from references 4 and 5. 


Description of Alternators Used for Program Evaluation 

Brayton cycle alternator . - The Brayton cycle alternator is rated 12 kilowatts at 
0. 8 power factor (lagging), 120/208 volts, 400 hertz, and 12 000 rpm. It is designed to 
be cooled with oil which has a temperature of 93° C. 

The stator laminations are 0. 007-inch electrical sheet steel and the yoke is made of 
ingot iron. Both the armature and field winding conductors are copper. The armature 
conductors are stranded and laid flat in the slot to minimize eddy-current losses. 

The rotor is made from AlSl 4620 steel and has laminated pole tips of 0. 014-inch 
electrical sheet steel. The laminated pole tips are electron-beam welded to the rotor 
and were used to minimize pole-face losses. In addition, zirconium copper damper bars 
were installed in the pole tips to equalize the terminal voltage during unbalanced loading 


r Damper bar-slot opening 


/ 



Figure 10. - Brayton cycle alternator rotor. 


C-68-3767 
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conditions. A photograph of the rotor is shown in figure 10. Note that some of the rotor 
material between the poles has been removed to reduce the leakage flux between poles 
from the rotor to the stator. Complete details of the alternator design are given in the 
sample output (pp. 19 to 26). 

SNAP- 8 alternators . - The two SNAP- 8 alternators are rated 60 kilowatts at 0. 75- 
power factor (lagging), 120/208 volts, 400 hertz, and 12 000 rpm. They are designed to 
be cooled with a polyphenyl ether oil, which has a temperature of 99° C 

A comparison of the magnetic materials used in the preprototsrpe and prototype 
SNAP- 8 alternators is shown in the following table: 



Preprototype 

Prototype 

Stator laminations (0. 014 in. ) 

AISI M- 19 

AISI M-19 

Rotor 

AISI 4130 

AISI 4620 

Yoke 

Ingot iron 

AISI 1020 


The prototype alternator has a thicker yoke than the preprototype and also has some 
of the rotor material between the poles removed as in the Brayton cycle alternator. In 
addition, the prototype alternator had circumferential grooves machined in the pole face 
surfaces in an attempt to reduce pole-face losses. Test results indicated that there was 
no major difference in pole-face loss between the prototype and preprototype alternators. 


Comparison of Experimental and Calculated Results 

Open-circuit satur ation curves . - A comparison of the test data and calculated re- 
sults for the open-circuit saturation curves of the three alternators are shown in fig- 
ure 11. In the corcq)uter program, field currents are calculated for a range of terminal 
voltages. The minimum voltage is 80 percent of rated terminal voltage. The volt^e is 
then increased by varying steps (maximum of 10 percent of rated terminal voltage) until 
some part of the magnetic circuit saturates. Saturation occurs when a flux density in a 
part of the circuit exceeds the maximum flux density of the appropriate material as spe- 
cified in the material data deck. 

The maximum percent difference between the experimental and calculated field cur- 
rents for the three alternators is 7 percent over the range of voltages from minimum to 
maximum. At rated voltage, the maximum difference is 4 percent. 
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(a) 400-Hertz Brayton cycle alternator. 




Field current, A 

(b) SNAP-8 alternator (preprototype). (c) SNAP-8 alternator (prototype). 

Figure 11. - Alternator open-circuit saturation curve. 


Field currents under load . - In table III, field currents are compared at rated volt- 
age and power factor for various alternator loads. At rated load conditions, the maxi- 
mum percent difference between the test and calculated field currents for any of the al- 
ternators if 5 percent. 

Losses and efficiency . - Before discussing losses and efficiency, the test and calcu- 
lated values for the field and armature winding resistances will be compared. This is of 
interest because, in determining copper losses, it is important that the winding resist- 
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TABLE in, - COMPARISON OF EXPERIMENTAL AND CALCULATED FIELD 


CURRENTS AT RATED VOLTAGE AND POWER FACTOR 


Alternator 

Load 

Field current, A 

Percent dif 


(a) 

Test'’ 

Calculated 

ference 

Brayton cycle (400 Hz) 

7. 5 kVA at 0. 8 power factor 

°4.1 

4. 1 

0 


15. 0 kVA at 0. 8 power factor 

5.7 

5.6 

1.9 


22, 5 kVA at 0. 8 power factor 

7.7 

7. 5 

2.6 


30, 0 kVA at 0. 8 power factor 

9.8 

9. 5 

3. 1 

SNAP- 8 (preprototype) 

60 kVA at 0. 75 power factor 

"^23. 1 

23.8 

3.0 

SNAP- 8 (prototype) 

60 kVA at 0. 75 power factor 

‘'l9. 1 

18. 1 

5.4 


^All power factors are lagging. 
^For separate excitation. 

^Test values from ref. 3, 
^Test values from ref. 4, 


ances be computed accurately from the conductor size and physical dimensions of the 
coil. 

A comparison of test and calculated winding resistances for the three alternators at 
25° C is given in table IV. All the corresponding test and calculated resistances agree 
to within 5 percent except for the calculated SNAP- 8 preprototype armature resistance 
which is low by 10 percent. 

The reason for this larger error is probably as follows. When the cross-sectional 


TABLE IV. - CALCULATED AND EXPERIMENTAL 
WINDING RESISTANCE AT 25° C 


Alternator 

Winding 

Resistance, ohms 

Percent dif- 
ference 



Test 

Calculated 

Brayton cycle (400 Hz) 

Armature 

Field 

®-0. 0382 
3. 27 

0.0389 
3. 30 

1.8 

.9 

SNAP- 8 (preprototype) 

Armature 

Field 

’’O. 0063 
1.46 

0. 0057 

1. 53 

10.0 

4.7 

SNAP- 8 (prototype) 

Armature 

Field 

'’O. 0057 
1.48 

0. 0056 

1. 53 

1.8 

3.3 


^Test values from ref. 2. 
'^Test values from ref. 4. 
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area of a rectangular conductor is determined in the program, the radius of the rounded 
corner is calculated per ASTM B48-55. The armature conductor of the SNAP- 8 prepro- 
totype alternator appears to have a larger corner radius than that used in the program. 
Hence, the computed conductor cross-sectional area is greater than the actual value. 
This results in a lower calculated than actual value for this particular resistance. 

Test and calculated values of the losses and electromagnetic efficiency at rated load 
and power factor for each of the three alternators is given in table V. For the test data, 
the method of separation of losses as given in reference 14 was used. For comparison 
of loss data, the following experimental losses are used: field and armature conductor 
losses and open-circuit core, and stray load losses. Since these are not the losses spe- 


TABLE V. - COMPARISON OF EXPERIMENTAL AND CALCULATED LOSSES 
AND EFFICIENCY AT RATED LOAD 


Alternator 

Load 

Loss or efficiency 

Test data 

Calculated 

Percent dif- 


(a) 

being compared 


data 

ference 

Brayton cycle (400 Hz) 

15 kVA at 0. 8 

Armature conductor, W 

‘=277 

272 

1.8 


power factor 

Field conductor, W 

135 

139 

2.9 



Open-circuit core, W 

320 

286 

11. 2 



Additional load,^ W 

270 

260 

3.8 



Total loss, W 

1002 

957 

4.6 



Efficiency, percent 

92. 3 

92.6 

. 3 

SNAP- 8 {preprototype) 

60 kVA at 0. 75 

Armature conductor, W 

‘^1470 

1344 

9.0 


power factor 

Field conductor, W 

1210 

1337 

10.0 



Open- circuit core, W 

1250 

1335 

6.6 



Additional load,^ W 

2500 

1995 

22.4 



Total loss, W 

6430 

6011 

6.7 



Efficiency, percent 

90. 3 

90.9 

.7 

SNAP- 8 (prototype) 

60 kVA at 0. 75 

Armature conductor, W 

®1320 

1323 

0.2 


power factor 

Field conductor, W 

800 

744 

7. 2 



Open circuit core, W 

1250 

1314 

5.0 



Additional load,^ W 

2000 

1883 

6.0 



Total loss, W 

5370 

5264 

2.0 



Efficiency, percent 

91.8 

91.9 

. 1 


^All power factors lagging. 

^Stray load loss for test data. Total of stator copper eddy, miscellaneous load and additional pole 
face, damper, and stator tooth and core due to load for calculated data. 


^Test values from ref. 3. 
^Test values from ref. 5. 
®Test values from ref. 4. 
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cifically calculated in the program, to make a comparison with the test data, some of the 
computed losses had to be added together. A table that shows the calculated losses cor- 
responding to the experimental values of the open-circuit core and stray load losses 
follows. 


Experimental loss 

Corresponding calculated losses that 
are added together 

Open-circuit core 

No-load pole face 
No-load stator tooth 
No-load stator core 
No-load damper 

Stray load 

Armature conductor eddy 
Miscellaneous load 
Additional pole factor, stator tooth, 
stator core, and damper due to load 


The maximum difference between the test and calculated values of electromagnetic 
efficiency for any one of the three alternators was 0. 7 percent. Agreement between the 
test and calculated data for the specific losses was not as good, ranging up to a maximum 
difference of 22 percent. Conductor losses can be in error due both to inaccuracies in 
the resistance computation and in the estimated operating temperature of the windings. 
The accuracy of the pole-face, tooth, and core loss calculations, all of which are highly 
empirical, affect the comparisons for the open-circuit core losses and for the additional 
losses due to load. 

Experimental and calculated values of electromagnetic efficiencies for the Brayton 




(a) 400-Hertz Brayton cycle alternator. 

Rgure 12. - Electromagnetic efficiency. 


(b) SNAP-8 alternator (preprototype). 
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cycle alternator over a range of loads from 25 to 125 percent of rated load are given in 
figure 12(a). Figure 12(b) shows a similar comparison for the SNAP- 8 preprototype al- 
ternator up to rated load. Maximum difference in data for the Brayton cycle alternator 
is 0. 8 percent. For the SNAP- 8 alternator, the maximum difference in test and calcu- 
lated efficiencies is 2 percent which occurs at 25 percent of rated load. From 50 percent 
to rated load, the maximum difference is 1. 0 percent. The difference at lower loads is 
not due to a large error in any one particular calculated loss. Rather, it is caused by an 
accumulation of small errors in several of the computed losses. 

Reactances . - A limited evaluation of the accuracy of the alternator reactance cal- 
culations was made. The direct-axis synchronous, and direct-axis transient reactances 
of the alternators were the only ones for which both experimental and calculated values 
were available. A comparison for these reactances is given in table VI. Except for the 


TABLE VI. - EXPERIMENTAL AND CALCULATED ALTERNATOR REACTANCES 


Alternator 

Reactance 

Test value, 

Calculated value, 

Percent dif- 



per unit 

per unit 

ference 

Brayton cycle (400 Hz) 

Direct-axis synchronous 

^1. 19 

1.21 

1.7 


Direct-axis transient 

.475 

. 563 

17.0 

SNAP- 8 (preprototype) 

Direct-axis synchronous 

*’l.40 

1. 57 

11.4 

SNAP- 8 (prototype) 

Direct-axis synchronous 

^1. 52 

1. 56 

2. 6 


Direct-axis transient 

.60 

. 656 

8.9 


®'Test values from ref. 3. 
^Test values from ref. 5. 
*^Test values from ref. 4. 


test data of the transient reactances, all values of reactances are for unsaturated condi- 
tions. 

The maximum difference between the experimental and calculated data for the 
direct-axis synchronous reactance is 11 percent. This is for the SNAP- 8 preprototype 
alternator. For the other two alternators, agreement is much better, being within 3 per- 
cent. For the direct-axis transient reactance, the calculated values exceed the corre- 
sponding test values by as much as 17 percent. This is probably due mainly to neglecting 
the effects of saturation on the calculated value. 


CONCLUDING REMARKS 


This report presents a digital computer program which calculates the electrical per- 
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formance characteristics of a homopolar inductor alternator from design data. A com- 
parison was made between the test results and calculated data for the 400-hertz Brayton 
cycle and SNAP- 8 alternators. The following observations were made. 

1. For the open-circuit saturation curves, the maximum difference between the test 
and calculated values of field currents was 7 percent. 

2. At rated load and power factor, the test and calculated field currents agreed to 
within 5 percent. 

3. The calculated efficiencies of the alternators at rated load and power factor were 
in agreement with the test results by a maximum difference of 0. 7 percent. 

4. For a range of alternator loads from 25 to 125 percent of rated load, test and cal- 
culated efficiencies agreed to within 2 percent. 

The program accuracy, as summarized above, is sufficient to allow using the pro- 
gram in practical applications such as parametric system studies and for specific alter- 
nator designs. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, October 28, 1968, 

120-27-03-42-22. 
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APPENDIX A 


COMPLETE FORTRAN LISTING AND FLOW CHARTS OF HOMOPOLAR 
INDUCTOR ALTERNATOR COMPUTER PROGRAM 

The complete FORTRAN listings of the main program and the three subroutines, 
which together constitute the homopolar inductor alternator computer program, are con- 
tained herein. The main program is INDCT, and the three subroutines are, in the order 
given, SINDUC, MAGNET, and OUTPUT. Each program listing, except that for OUTPUT, 
is followed by its flow chart. The organization of OUTPUT is self-evident since it con- 
sist largely of WRITE and FORMAT statements. 


INDCT 

INCCT 


C 

C 


C 


1 

2 

3 

C 

C 

C 

4 

5 

6 

7 

8 
C 


CONNCN 4,4^ ,AB,ACt ACRf AGfAI, ALh.ALPhAEt ALPHAR,ALPHAS,ALY,ALYC»ALYR A 
I ,AP,AS, ASH, ATM, AY, AYC, AYR, B,B 1,82, B3.8 CLL, BCCIL.bg .BK, BN, 80, BP, BPL A 
2.8S,BSHL,BTLL,BV.BYCLL,C ,C1 , CC , CCR ,CE ,CF ,CK ,CL , CM , CP ,CQ , CW , D, D 1 , DC A 
3CIL ,0D,DF,C I ,DISH,CISH l,CR,DSH,DU.DW,DWl,DYC,EC.EDD,EE.EL,EP,Eh.F, A 
AFCL ,FE,FFLL,FGL,FGML,FF,FK1,FPL,FC ,FS,FSHL,FSHLP,FTL,FYCL,FYL,FYRL A 

5 , G,GA,GC,GE,GP ,GXX,H,FC,hD,HM,Ha,HP,HPl.hS.HT,HV,HW,HX,HY, IBN, IPN, A 
6 IPX. IQQ,I2Z, JA,KSAT,LTR,LTR1 , L TS , P5 , P6 , P7 , PB A , PC , PCC I L , PE , PF , PHL , P A 
7HU,PI,PL,PF,PML,PN,PT ,PX ,QN , CQ . RC , RD . RE . RF ,RG 1 . RK , RK 1 , RPM , RR , RS , RT A 

8, RY,$,Se,SC,SC,SF,SH, S I, SI6MA,SK,SN, SN1,SS,STATET,T1,T11 ,T2,T22,T3 A 

9, T33,TB,TC,TF,TG.TS. TST.TT.TY.TYE.TYPY.TYR.VA. VR. WC.WF.WI .WL.WO WR A 

SCTCR ,WTCTAL,V<YOKE,XA, XB,XD,XF,XL,XQ,XR,XU,YY,Z,ZG,ZZ,ZZZ A 

A 

INTEGER TYPY,ZZ A 

A 

CIMENSICN QVLNUO), QVLL(IO), GFCUR(IO), QAGATllO), QTATUO). QPAT A 
1(10, CCAT(IO), CTFAT(IO), QSAT(IC), QYAT(IO), QPD(IO), gCO(lO), Q A 
27FC(I0), CSC(IO), CYO(IO), GPERV(IO), AK90), GX ( 5 ) , YA(5), ED{5), A 
3 FGX(5), PR(5), FI(?>). PS(5), G(5), DL(5), PP(t>), EX{5). ST(5), WA A 
A(5), STRAY(5), FFL(5), BSHLL(5), BCL(5), BTL(5). BPLL(5). BYCL(5), A 
5 FGLLv5). PFLL(5), P(5), £(5), PZ(5), SP(5), WGL(5), CDD(5), EF(5) A 

6, AKVA(5), EZ(IO), TTA(IC), TTB(IO), RRA(IO), RRB(IO), SMAT(6), RM A 


7AT(6), YMAT(6) A 

A 

READ (5,1) SPAT A 

READ (5,2) (AI(I),I=l,29) A 

READ (5.1) RMAT A 

REAC (5,2) (AI(I),I^31,59) A 

REAC (5,1) YPAT A 

REAC (5,2) (Aid), 1=61, 69) A 

FORMAT (6A6) A 

FORMAT (8F10.1) A 

CALL SINDUC A 

CALL OUTPUT A 

A 

CCMFUTE TCCTH WIDTH AT 1/3 DISTANCE FROM NARROWEST SECTION A 

A 

IF (ZZ-3) A, 5,6 A 

SM=TT-8S A 

GC TO 8 A 

SM= (3.1A16«(Cd2.*FS)/CG )-B3 A 

GO TO 8 A 

IF (ZZ-A) 5, 7, A A 

SM=TT-.9A*eS A 

CCNTINUE A 


A 


1 

2 

3 

A 

5 

6 

7 

8 
9 

10 

11 

12 

13 

lA 

15 

16 

17 

18 

19 

20 
21 
22 
23 
2A 

25 

26 

27 

28 

29 

30 

31 

32 

33 
3A 

35 

36 

37 

38 

39 
AO 
A1 
A2 
A3 
AA 
A5 


INDCT 


35 



c 

ARE/JS AND LENGTHS FOR MAGNETIC CALCULATIONS 

A 

46 

c 



A 

47 


AP=8P*PL*RK 


A 

46 


ACR=(DU-2.*HC)*3.1416«PE*SS/PX 


A 

49 


FATI=FGPL 


A 

50 


ATH=GQ*SS*SN*PE/PX 


A 

51 


ASH=<DSH*«2-CISH«»2)*,7854 


A 

52 


AY=TY*(CU+TY)*3.1416 


A 

53 


IF (TYPY-2) 10,11,9 


A 

54 

9 

ALY==1.334*CL 


A 

55 


GC TO 12 


A 

56 

10 

AYR = 0 


A 

57 


AYC = 0 


A 

58 


ALY=BCCIL+.667*CL 


A 

59 


ALYR=0 


A 

60 


ALYC=0 


A 

61 


GO TO 13 


A 

62 

11 

ALY-.667*CL 


A 

63 

12 

AYC^3.1416* (CYC+TYE)*TYE 


A 

64 


AYR=TYR*(CU+2.*TY)«3. 1416 


A 

65 


ALYC=BCCIL 


A 

66 


ALYR=*OYC-DL 


A 

67 

13 

CCNTINLE 


A 

68 

C 



A 

69 

c 

NO-LOAD, RATED VOLTAGE MAGNETIZATION CHARACTERISTICS 

A 

70 

c 



A 

71 


22Z=PX»GE/( .C0319*GA«PE) 


A 

72 


KSAT=10 


A 

73 


GXX=1. 


A 

74 


ECC=1. 


A 

75 


FH^eG*GE/0.C0319 


A 

76 


FGML=0. 


A 

77 


CALL MAGNET 


A 

78 


Jsl 


A 

79 


FGLL(J)=FGL 


A 

80 


PNLL(J)=PML 


A 

81 


BPLL(J)=BPL 


A 

82 


eTL(J)=BTLL 


A 

83 


BSHLLl J )=eSHL 


A 

84 


BCL( J)=BCLL 


A 

85 


BYCL(J)=^BYCLL 


A 

86 


FFLU)=FFLL 


A 

87 

c 



A 

88 

c 

SHORT CIRCUIT RATIO AND SHORT CIRCUIT AMPERE-TURNS CALCS 

A 

89 

c 



A 

90 


FSC=XA»FH*C.02 


A 

91 


SCR=FFLL/FSC 


A 

92 


WRITE (6tl4) FSC.SCR 


A 

93 


FORMAT (1HL,9X,27H SHCRT CIRCUIT 

AMPERE-TURNS, F16*3/10X,20H SHORT 

A 

94 


ICIRCUIT RATICfF23.3) 


A 

95 


WRITE (6,15) SMAT 


A 

96 

15 

FORMAT (1HL,18H STATOR MATERIAL - 

,1H ,6A6) 

A 

97 


WRITE (6,16) RMAT 


A 

98 

16 

FORMAT (1HL,18H RCTOR MATERIAL -- 

,1H ,6A6) 

A 

99 


WRITE (6,17) YMAT 


A 

100 

17 

FORMAT (1HL,18H YOKE MATERIAL 

,1H ,6A6) 

A 

101 
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FYOKE=FYL+FYCL*FYRL A 102 

KRITE 16 , 18 ) TG,Fe,BG,BPL, BILL, BCLL.BSHL.BYCLL.FGL, FPL, FTL.FCL.FSH A 103 
1LP,FSHL,FYCKE,FFLL A lOA 

18 FORPAT (1H1,30H MAGNETIZATION CHARACTER I ST IC S/ 5X ,25H (NO LOAD, RAT A 105 

lED V0LTAGE)//10X,18H TOTAL USEFUL FLUX , F 12.2 , lOH K ILOL INES/lOX , 1 7H A 106 
2 USEFUL FLLX/P0LE,F13.2//10X,15H FLUX DENSIT IES/13X, 19H AIRGAP (IN A 107 
3CL. PHL),F8.2,9H KL/ SC- IK/ 13X , 5H P0LE,F22.2/ I3X,6H TOOTH ,F21 . 2/13X A 108 
A,5H C0RE,F22.2/13X,19H SHAFT (UNDER FLO. ) , F8.2/ 1 3X , 17H YCKE (OVER A 109 
5FLD.),F10.2//10X,13H AMPERE-TURNS/ 13X ,7H A IRGAP , F20.2 , 1 IH PER STAT A 110 
60R/13X,5H P0LE,F22.2/13X,6H TCOTH , F2 1 .2/13X, 5H CORE , F22 .2/ 13X, 19H A 111 
7SHAFT (UNDER PCLE > ,F8 . 2/ /I 3X 19H SHAFT (UNDER FLD. ) ,F8.2/13X, 5H YOK A 112 
8E,F22.2//13X,6H T0TAL,F21.2) A 113 

IF (KSAT.EC.O) GO TO 19 A 114 

GO TO 20 A 115 

19 WRITE (6,76) A 116 

GO TO 3 A 117 

C A 118 

C HCT AND COLD CAMPER BAR LOSS CALCULATIONS A 119 

C A 120 

20 IF (BN) 21,21,22 A 121 

21 WO=C.O A 122 

KU=C.O A 123 

GO TO 44 A 124 

22 AA=hO/GE A 125 

VT=C A 126 

IF (AA) 23,26,23 A 127 

23 IF (AA-C.65) 24,26,25 A 128 

24 VT=AL0G(10.«AA)*(-0.242)+0.59 A 129 

GO TO 26 A 130 

25 VT=»0.327-(AA*0.266) A 131 

26 CONTINUE A 132 

FS1=2.0»QN*PN*F A 133 

FS2=2.0*FS1 A 134 

M=0 A 135 

RM=RE*( 1.CT4ALPHAE*(T33-2C. ) ) A 136 

GO TO 28 A 137 

27 RM=RE*(1.04ALPHAE«(T3-20.) ) A 138 

28 AA= (FS1/RM)«*0.5*CC*0.32 A 139 

AB=(FS2/RM)««0.5*DD*0.32 A 140 

IF (AA-2.5) 29,29,30 A 141 

29 V1=1.0-0.15»AA40.3«AA«AA A 142 

GC TO 31 A 143 

30 V1=AA A 144 

31 IF (AB-2.5) 32,32,33 A 145 

32 V2=1.0-C.15*AB4C.3«AB*AB A 146 

GC TO 34 A 147 

33 V2=AB A 148 

34 IF (H.EC.O.) GO TO 35 A 149 

IF (H.EO.E) GO TO 35 A 15C 

VC=H/(3.0«E»V1) A 151 

60 TO 36 A 152 

35 VC=C.75/V1 A 153 

36 VS=HD/NC+VT4VC A 154 

VG=TB/(CC*GC) A 155 

01=1.0-(1.0/(((BO*0.5/GC)*«2.0+1.0)**0.5)) A 156 

CZ=EO/TS A 157 
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C2=^1.05»SINCC2«2.844) 

A 

158 


IF (QZ-0.37J 37,37,38 

A 

159 

37 

C3=C.46 

A 

160 


GC TO 39 

A 

161 

38 

C3=C.23*SIN(I0.46*QZ-2,U+0.23 

A 

162 

39 

C4 = SIN(6.283«TB/TS-l.57l )*1.C 

A 

163 


C5=SINtl2.566*TB/TS-1.571)*l.O 

A 

16A 


IF (H) 41,40,41 

A 

165 

AO 

4B=C.785»DD*C0 

A 

166 


GC TO 42 

A 

167 

A1 

»8=F*B 

A 

168 

42 

W2=PX*BN*SB*RM»1.246/ (iB*1000. ) 

A 

169 


W3=(Q2/ (2.0*VS+< VG/Q4 ) ) ) «»2.0»V1 

A 

170 


W5=CQ3/(2.C«VS+(VG/Q5) n««2.C*V2 

A 

171 


KC= (TS*BG»C1»CC >»*2.0*V»2*(W3+W5) 

A 

172 


F=F + 1 

A 

173 


IF (M-1) 44,43,44 

A 

17A 

A3 

V)U = hO 

A 

175 


GC TO 27 

A 

176 

AA 

CONTINUE 

A 

177 

C 


A 

178 

C 

PCLE-F4CE Less CALCUL4TICN 

A 

179 

C 


A 

180 


GT=eO/GC 

A 

181 


<jS=1.75/(GT«*1.35)+0.8 

A 

182 


GF=flA*PI*SC/tC*FH) 

A 

183 


C2=BG*»2.5*0.CC0061 

A 

18A 


D3=(0.0167*C0*RPM>**1.65»0.OC00l5147 

A 

185 


IF (TS-0.9) 45,45,46 

A 

186 

A5 

C4=TS*»1. 285*0. 81 

A 

187 


GC TO 49 

A 

188 

A6 

IF (TS-2.0) 47,47,48 

A 

189 

A7 

C4=TS**1. 145*0.79 

A 

190 


GC TO 49 

A 

191 

A8 

C4=TS**0. 79*0.92 

A 

192 

A9 

D7=B0/GC 

A 

193 


IF (07-1.7) 50,50,51 

A 

19A 

50 

C5=C7**2.31*C.3 

A 

195 


GO TO 56 

A 

196 

51 

IF (07-3. C) 52,52,53 

A 

197 

52 

05*07**2.0*0.35 

A 

198 


GC TO 56 

A 

199 

53 

IF (07-5.0) 54,54,55 

A 

200 

5A 

05=07**1.4*0.625 

A 

201 


GC TO 56 

A 

202 

55 

05=07**0.965*1.38 

A 

203 

56 

C6=10.0**(0.932*C1-1.6C6) 

A 

20A 


V.N = 01*C2*D3*D4*C5*C6*GA 

A 

205 

C 


A 

206 

C 

CALCULATE NC-LOAC , RATED VOLTAGE TOOTH AND CCRE LOSS 

A 

207 

C 


A 

208 


hT=(SM)«CC*SS*HS*0.849*( PTL( 1 )/BK)**2.0*WL 

A 

209 


hC=(DU-HC)*2.67*HC*SS« (BCL( 1 )/BK)**2.0*WL 

A 

210 
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c A 211 

C ARRANGING LCAO POINTS IN ORDER A 212 

C A 213 

DC 58 J=l,^ A 214 

IA=5-J A 215 

DC 58 1 = 1, lA A 216 

IF (G( n.GT.G(I-»in GC TC 57 A 217 

GO TO 58 A 218 

57 POL=GlI) A 219 

G(I)=G(I+1) A 220 

Gn*l) = POL A 221 

58 CONTINUE A 222 

Gm=0. A 223 

PN=5 A 224 

CO 59 1=2,5 A 225 

IF IG(I).GE.1.0.ANC.G(I-1).LT. 0.999) MM=I A 226 

YA(I) = 1C0./Gm A 227 

59 CONIINUE A 228 

C A 229 

C CALCULATE GENERATOR LOAD CHARACTERISTICS A 230 

C A 231 

AN=ARCCS(PF) A 232 

DO 60 J=2f5 A 233 

AA=4TAN( (Xe/YA( J )+SIN (4N))/PF) A 23A 

EB=4A-AN A 235 

EO(J) = XA*SIN(AA)/YAC J )+CCS(8B) A 236 

FGX(J)=FATI*lCO./YA( J) A 237 

GX( J)=( (EC(J)-C0*93*XC*SIN(AA)/YA( J) ) ))*CK A 238 

TTB(J)=0. A 239 

TTA(J)=C. A 240 

PRA(J)=C. A 241 

FRfi(J)=0* A 242 

E2(J)=0. A 243 

STRAY(J)=0 A 244 

PRLL(J)=0 A 245 

FFL(J)=C A 246 

BSHLL(J)=C A 247 

BCL(J)=0 A 248 

BTL(J)=C A 249 

BPLL(J)=0 A 250 

EYCL(J)=0 A 251 

FI(J)=0 A 252 

CCC(J)=0 A 253 

EF(J)=0 A 254 

PR(J)=0 A 255 

ST(J)=0 A 256 

V(CL{J>=G A 257 

FP(J)=0 A 258 

CL(J)=0 A 259 

PS(J)=0 A 260 

EXCJ)=0 A 261 

SP(J)=0 A 262 

AKV4(J)=0 A 263 

^4(J)=0 A 264 

P(J)=0 A 265 

PZCJ)=0 A 266 

E(J)=0 A 267 

60 FGLL(J)=0 A 268 

J*2 A 269 
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II 



KSAT^^IO 

A 270 

61 

GXX=6X( J) 

A 271 


ECD=ED( J) 

A 272 


FGML=FGX( J) 

A 273 


CALL HA6NET 

A 274 


FGLL(J)=FGL 

A 275 


PMLL(J)=P^»L 

A 276 


BPLL(J)=BPL 

A 277 


eTL(J)=BTLL 

A 278 


BSHLL( J)=BSHL 

A 279 


BCL( J)*BCLL 

A 280 


BYCL(J)=BYCLL 

A 281 


FFL< J)=FFLL 

A 282 


IF (KSAT.EQ.O) GO TO 62 

A 283 


IF (J.EC.5) GO TO 62 

A 284 


J*J + 1 

A 285 


GC TO 61 

A 286 

62 

JA = J 

A 287 


IF (KSAT.EC.O) JA*JA-1 

A 266 


FI (PM)*FFL(KP)/PT 

A 269 


ViVf^WU 

A 290 


VV=3**PI*EP*PF 

A 291 


P=1 

A 292 

63 

CONTINUE 

A 293 

C 


A 294 

C 

EDDY FACTOR CALCULATIONS 

A 295 

C 


A 296 


UA = 6(M) 

A 297 


TTA(M) = (T1--T11)*UA*UA4T11 

A 298 


RB=(1,0E-6)*RS*( 1.0+ALPHAS*(TTA(M)-20.) ) 

A 299 


IF (SH) 64,64»65 

A 300 

64 

EZ(P)=1. 

A 301 


GC TO 66 

A 302 

65 

AA=0.584+(SN*SN-1.0)*C*0625*(S0*CL/(SH*HM/2.) )**2 

A 303 


AB=(SH*SC»F*AC/(BS*RB*10C0000.0) )**2,0 

A 304 


ET= A A* A 8*0.00335+1.0 

A 305 


EB=ET-0.00168*AB 

A 306 


EZ(N)=(ET+EB)*0.5 

A 307 

C 


A 308 

C 

LOSSES AND EFFICIENCY UNDER LOAD 

A 309 

C 


A 310 

66 

FI (N) = FFL(N)/PT 

A 311 


CDD(M)=FI (P)/AS 

A 312 


TT8(M)=nT2-T22)/CFI (PN)-FK 1) )*»2)*(FI(M)-FI(1) )**2+T22 

A 313 


RRB(M)= (1.0E-6)*RR*( 1.0+ALPHAR*(TTB(M)-20. ) )*ZG 

A 314 


PR (P)=FnN)*FHNMRRB{M) 

A 315 


EF(P)*FI(M)*RRB(M) 

A 316 


RRA(M)-R6«RY 

A 317 


PS(P)s(3.«(PI*UA)**2)*RRA(H) 

A 316 


MOL(M)»WQ*(BCL(H)/BCL(n )*«2 

A 319 


ST(P)=WT*(BTL(N)/BTL( 1) ) **2 

A 320 


KA(P)*VV*UA/1000. 

A 321 


AKVA(M)=WA(M)/PF 

A 322 


STRAY(M)=AKVA(M)*10.0 

A 323 


GM=(GF«UA)«*2. 0+1.0 

A 324 


CL(R)-GP*kW 

A 325 
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o u u 


67 


68 


69 


PP(^)=G^'*WN 

EX (P)* (EZ(P )-U0 )*PS(^ )*2*C*CL/HM 

SP(P)=PP(M)+OL(M)+PR(K )+PS(M)+EX(M)+ST(M)+WF+WQL(M)+$TRAY(M) 
P(M)=(SP(M)/1C00. )+WA(P) 

PZ(^*)=(SP(P)/P(H) )*.l 
E(f*) = lCO.C-PZ(f') 

IF (M.EC*1) UW=WO 

IF (M.LE.JA) GC TC 63 

WRITE (6,67) PF, (G(I) ,1=1,5) , (PMLLd ),I = l,5) , (FGLL ( I ) , 1 = I • 5 ) , (BPLL 
l(I),l=l,5),(BTL(I),I=lt5), (BSHLLU ), 1=1,5) , ( BCL ( I ) • I=l ,5 ) • ( BYCL (I ) 
2,1 = 1,51 ,(FFL( 1 1, 1=1, 5), (FI (I ), 1 = 1,5 ),(CDO( I), 1=1,5 ),(EF( I), 1 = 1, 5) 
FCRPAT (lHl,26XA7H/iLTERN/JTCR LCAO CHARACTERISTICS (RATED VOLTAGE, F 

15.2, 14H POWER FACTOR ) /27X66H 

2 ////7X12HPERCENT LOAD, IIX , 2PF17.0,4F 19 

3.0//7X18HLEAKAGE FLUX ( PPL ) , 5X , OP5F 19 .2/7X2CHA IR-GAP AMPERE TURNS, 
43X5F19.2//7X25HFLUX DENSITIES (KL/SQ-IN ) /ICXAHPOLE , 16X5F19,2/10X5H 
5TEETH15X5F19.2/10X18HSHAFT (UNDER FLD. ) 2X5F 19, 2/ 10X4HCCRE 16X5F 19. 2 
6/10X16HYOKE (OVER FLD . ) 4X5F 19. 2//7X18HTOTAL AMPERE TURNS5X5F19.2/7 
7X20hFIEL0 CURRENT (AMPS ) 3X5F19.2/7X21HCURRENT DENS. ( F lELD ) 2X5F19. 
82/7X11HFIELO VOLTS 12X 5 Fl 9. 2 ) 

WRITE (6,68) (TTB(I),I=1,5), (TTA(I),I=1,5) ,(RRB(I),I=1,5),IRRA(I), 
11=1,5) , (EZ( I) , 1=1,5) , (PR(I ), 1 = 1,5) ,WF,WF,WF,WF,WF, (ST (I) ,1 = 1,5), (W 
2CL(I),I = 1,5),(PP(I ),I = l,5),(CLm.I = l,5).(PS(l ),I=1,5),(EX(I),I = 1, 
35) , (STRAY ( I ) ,1 = 1,5) , (SP( I ) , I = 1 , 5 ) , ( AKVA ( I ),I = 1,5),(WA(I),I = 1,5),(P 
4(1), 1=1, 5), (PZ (I), 1=1, 5), (Ed), 1 = 1, 5) 

FORMAT ( 1HK,6X20HTEMPERATURES ( DEG .C ) / IOX5HF lELD 15X5F 19 . 2/ 10X8HARM 
1ATURE12X5F19.2//7X18HPES ISTANCES ( OHMS ) / 10X5HF lELO 15X5F 19 .2/ 10X8HA 
2RMATURE12X5F19.^//7X11FEC0Y FACTOR 12X5F 19. 2//7X25HALTERNA TOR LOSSE 
3S (WATTS)/10X5HFIELD15X5F19.2/10X7HWIN0AGE13X5F19.2/10X12HSTATOR T 
4CCTHaX5F19.2/10XllHSTATGR C0RE9X5F 19 .2/ 10X9HPOLE FACE 1 1X5F 19 . 2/lOX 
56HCAMPER14X5F19.2/10X13HSTATCR COPPER 7X 5F 19 • 2/ 10X4HEDD Y 16X5F 19 .2/ 1 
60X1CHMISC. LCAC10X5F 19.2/1CX5HTUTAL15X5F 19 .2// 7X23HALTERNATOR OUTP 
7LT (KVA )5F19.2/7X22HALTERNATOR OUTPUT ( KW ) 1X5F19. 2/7X2 IHALTERNATOR 
8 INPUT (KW)2X5FI9.2/7X14FPERCENT LGSSES9X5F 19. 2/ 7X 18HPERCENT EFFIC 
9IENCY5X5F19.2) 

CALCULATE NO-LOAD SATURATION DATA 

CO 69 J=l,10 
CPERV( J)=C 
CVLL( J )=0 
CVLN(J)=0 
CFCUR( J)=0 
CTAT(J)=0 
CAGAT( J)=0 
CPAT( J)=0 
CCAT(J)=0 
CTHAT( J ) = 0 
GSAT(J)=0 
CVAT(J)=0 
CPD(J)=0 
CCD(J)=0 
GTHC(J)=0 
QSD( J)=0 
CYC ( J)=C 


A 326 
A 327 
A 328 
A 329 
A 33C 
A 331 
A 332 
A 333 
A 334 
A 335 
A 336 
A 337 
A 338 
A 339 
A 340 
A 341 
A 342 
A 343 
A 344 
A 345 
A 346 
A 347 
A 348 
A 349 
A 350 
A 351 
A 352 
A 353 
A 354 
A 355 
A 356 
A 357 
A 358 
A 359 
A 360 
A 361 
A 362 
A 363 
A 364 
A 365 
A 366 
A 367 
A 368 
A 369 
A 370 
A 371 
A 372 
A 373 
A 374 
A 375 
A 376 
A 377 
A 378 
A 379 
A 380 
A 381 
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FGHL=0. A 382 

FG=FH A 383 

ECC^l* A 38A 

GXX=1. A 385 

IDELR=10 A 386 

R=.e A 387 

J=1 A 388 

KSAT=10 A 38S 

FH=FG*R A 39C 

FG=TG*R*CP/PX A 391 

CALL MAGNET A 392 

IF (KSAT.EC.O) GC TO 71 A 393 

CPERVr J)=1C0.«R A 39A 

CVLL(J)=EE*R A 395 

GVLN(J)=QVLL( J)/SQRT(3.) A 396 

CFCLR( JMFFLL/PT A 397 

CTAT(J)^FFLL A 398 

CAGAT(J)=FGL A 399 

CPAT(J)*FPL A 400 

CCAT(J)=FCL A 401 

CTFAT(J)=FTL A 402 

CSAT(J)^FSHL+2.*FSFLP A 403 

CYAT(J)=FYL+FYCL+FYRL A 404 

CPC(J)=BPL A 405 

CC0(J)=8CLL A 406 

CTHC(J)=BTLL A 407 

CSD(J)=BSFL A 408 

CYC(J)=BYCLL A 409 

IF (J^EG.IC) GO TO 74 A 410 

J=J+1 A 411 

R:=R + FLGAT ( I0ELR)/1C0. A 412 

GC TO 70 A 413 

B=R-FLCAT( ICELR)/100. A 414 

IF (IDELR.EC.IO) GO TC 72 A 415 

IF (I0ELR.EC.5) GC TO 73 A 416 

GC TO 74 A 417 

ICELR=5 A 418 

R = R + FLCATnDELR)/lCO. A 419 

KSAT=10 A 420 

GC TO 70 A 421 

ICELR=2 A 422 

R=R+FLCAT( IDELR)/1C0. A 423 

KSAT=10 A 424 

GC TO 70 A 425 

WRITE (6,75) (CPERV(K),K=1.10),(QVLN(K),K=1,10).(QVLL(K).K=1.10),( A 426 

1QFCUR(K),K=1,10), (CPD(K) ,K= 1 , 10 ) , ( QTFO ( K ) , K= 1 , 10), ( QSD ( K ) , K= 1 , 10 ) , A 427 

2(CCC(K) ,K=1,10) , (CYD(K ),K=1, 10 ) , ( C AGAT ( K ) . K= 1, 10 ) , ( GPAT ( K ) ,K= 1 , 10 ) A 428 

3, (CTHAT(K),K=1, 10) ,(QCAT(K),K=1,10) , ( QSAT ( K ) ,K= 1 , 10 ) , (GYAT(K) ,K=1, A 429 

410),(QTAT(K),K=1,10) A 430 

FORMAT (1H1,50X23HN0-LCAC SATURATION 0ATA/51X23H A 431 

1 //2X7HVCLTAGE/5X7FPERCENT6X10F11.2//5X12HLINE-NEUTRAL1X10F11* A 432 

22/5X9HLINE-LINE4X10F11.2//2X13FFIELD CURRENT3X lOF 1 1 • 2//2X20HFLUX 0 A 433 

3ENS. {KL/SC-IN)/5X4HP0LE9X10F11.2/5X5HT00TH8X10F11.2/5X5HSHAFT8X10F A 434 

411.2/5X4HCCRE9X10F11.2/5X4HYCKE9X10F11.2//2X12HAMPERE^TURNS/5X6HAI A 435 

5RGAP7X1OF11.2/5X4HPOLE9X1OF11.2/5X5HTOOTH8X10F11.2/5X4HCCRE9X1OF11 A 436 

6.2/5X5HSHAFT8X10FlU2/5X4HYOKE9X10F11.2//5X5hTOTAL8X10F11.2) A 437 

GC TO 3 A 438 

FORMAT (lHL10X,17hMACFINE SATURATED) A 439 

END A 440- 
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Calcuiate tooth wictth 
at 1/3 distance from 
narrowest section (SM) 


'WRITE: 

Rating 
slots 
airgap 

armature winding 

field winding 

stator 

rotor 

shaft 

damper bars 
yoke 
weights 
constants 
permeances 


8 >13 
Calculate areas and 
lengths for mag- 
netics calculations 


No-load rated-voltage 
magnetization characteristics 


Calculate short- 
circuit ratio fSCR) 
and ampere - turns 
IFSC) 


14-18 




KSATy 


0-0 


Input-output 


V 


Branch pointer decision 


□ 


Arithmetic statement 
Call statement (or 
subroutine 


Numbers In upper left indicate external 
formula number in FORTRAN program 


^0 
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SUBROUTINE SINCUC B I 

CCRRCN «(A«,/)8,AC>ACR.^G,AI.i)LH«ALPHAE,ALPHAR»ALPHAS,ALY,ALYC>ALYR B 2 

1. AP.ASt ASH. ATHf AY. AYC. AYR. 6. 81,62. B3.BCLL« BC01 L.bg tBK. BN »BO. BP. BPL B 3 

2, BS,BSHL,BTLL,BV,BYCLL ,C,Cl,CC,CCR,CE,CF,CK,CL.CM,CP,CQ,CVl,D,Dl,OC B A 

3CIL.OO.OF,CI,OISH.CISHl.CR.OSH,OU,DN,OWl.CYC,EC,EOO,EE.EL,EP.Eta.F. B 5 

AFCL.FE.FFLL.FGL.FGML.FF.FKl.FPL.FQ.FS.FSHL.FSHLP.FTL.FYCL.FYL.FYRL B 6 

5,G.GA,GC.6E.GP.GXX,H,FC.H0,HP.HC.HP,HP1,HS,HT.HV.HH.HX.HY,1BN. IPN, B 7 

eiPX.IQQ.IZZ, JA,KSAT,LTR,LTR1,LTS,P5,P6,P7,PBA,PC,PC0IL,PE,PF,PHL,P B 8 

7HW,PI,PL,PP,PNL,PN,PT,PX,QN,CC1,RC,RD,RE,RF,RG1,RK,RK1,RPP,RR,RS,RT B 9 

8, RY,S,SE,SC,SD,SF,SH,SI,SIGMA,SK,SN,SNl,SS,STATET,Tl,Tll,T2,T22,T3 B 10 

9, T33,TB,TC,TF,TG,TS, TST.TT.TY.TYE.TYPY.TYR.VA, VR.WC.WF.WI.WL.WO.WR B I I 


$CTOR,WTCTAL,WYOKE,XA,XB,XD,XF,XL,XQ,XR,XU,YY,Z,ZG,ZZ,ZZZ B 12 

C fl 13 

INTEGER TYPY.ZZ B lA 

C B 15 

BEAL LT.LTS.LTR.LTRl B 16 

C B 17 

DINENSICN CA18), DXf6), CY(8), DZ{8), AI(90>, G(5) B 18 

C B 19 


NAPELIST /RATING/ VA , E E , EP , F , RPM, I PX ,PF ,G/ST ATOR/D I .DU .CL , HV , BV , SF B 20 
1,LTS,WL,BK/SLCTS/ZZ,BC,B3,BS,H0,HX,HY,HS,HT, ICG/K INONG/RF , SC . YY, C , B 21 
2DW,SN,SNl,Chl,CE,SD,PEA,SK,Tl,RS,ALPHAS,Tll.TST/AIRGAP/GC,GP/C0NST B 22 
3/Cl, CP, EL,CN,CC,WF/R0TCR/RK,PL,HP,HP1,PE,BP,WR0T0R,LTR,LTR1,RK1,PH 8 23 
4h,PHL,Dl/CAf<P£R/WC,HD,CD,H,a,BN,SB,TB,T33,RE,ALPHAE,T3/ShAFT/0SH,D B 24 
SISH.DISFl, ALH/YCKE/TYPY, TY.TYE.TYR.OYC/FIELD/PCOIL.DCOIL.PT.RD.RT, B 25 


6T2,BC0IL,TF,T22,RR.ALPHAR B 26 

C B 27 

DATA DA, OX, OY.CZ/C. 05, C.C72.C. 125, 0.165, 0.225, 0.438, 0.688, 1.5, 0.00 B 28 

10 124. 0. 0002 1.0. 00021, C. 00084,2*0. 00 189, 2*0. 000124,2*0. 00084, 0.0018 B 29 

29.0. 00335.0.00754.0.0302.3*0.000124,2*0.00335,0.00754,0.0134,0.030 B 30 

32/ B 31 

C B 32 

C1=C B 33 

RS=C.694 B 34 

PR*C,694 B 35 

RE=0.694 B 36 

ALPFAS*C.CC393 B 37 

ALPHAR=0. 00393 B 38 

ALPFAE=0.CC393 B 39 

T33=20. B 40 

TF=25. B 41 

TST=25. B 42 

GP=C. B 43 

RK1*0. B 44 

PE=C. B 45 

PHV«=0. 8 46 

PL=C. 6 47 

PHL=0. B 48 

LTR1=0. B 49 

D1=C. B 50 


48 


SINDUC 



P8A=60, 


B 

51 

SN*1.0 


B 

52 

CYC=0# 


B 

53 

SH = 0. 


B 

54 

DVil = 0 


B 

55 

CC = C. 


B 

56 

Ch=0 


6 

57 

CP=0 


B 

58 

EL = C 


B 

59 

CP = C 


6 

60 

G( 1 1=0. 


6 

61 

G(2)*0.75 


B 

62 

G (3 l^l.CO 


8 

63 

G(4)=1.25 


B 

64 

G(5)=1.50 


B 

65 

CQ = C 


B 

66 

PP = C 


B 

67 

P5 = C 


B 

68 

P6=0 


8 

69 

P7 = C 


B 

70 

KC = C. 


B 

71 

UF = 0 


B 

72 

TY = 0 


6 

73 

TYE = 0 


B 

74 

TVR = 0 


B 

75 

CYC^O 


6 

76 

EP=C. 


B 

77 

EE = C. 


B 

78 

IPN = 3 


B 

79 

PN = 3. 


B 

80 

IPX = 0 


B 

81 

F = 0. 


B 

82 

RPP=0. 


B 

83 

BP = C. 


B 

84 

SF = C. 


B 

85 

RK = C. 


B 

86 

LTS=0. 


B 

87 

LTR=0. 


6 

88 

hROTOR=C. 


B 

89 

HV = C. 


B 

90 

BV = C. 


B 

91 

BCCIL=0* 


B 

92 

h=0. 


B 

93 

SK = C 


B 

94 

WRITE (6.1) 


B 

95 

1 FCRPAT ( IF 143X3 3F*»HGPCPCLAR 

INDUCTOR ALTERNATOR*#) 

B 

96 

REAC (5, RATING) 


B 

97 

REAC (5.STATCR) 


B 

98 

REAC (5.SLCTS) 


B 

99 

REAC (5,WINCNG) 


B 

100 

REAC (5,AIR6AP) 


B 

101 

REAC (5.CCNST) 


B 

102 

REAC (5,RCTCR) 


B 

103 

REAC (5, DAPPER) 


B 

104 

REAC (5. SHAFT) 


B 

105 

REAC (5,YCKE) 


B 

106 

REAC (5, FIELD) 


B 

107 


SINDUC 
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I 



o) o n n 


IF (EP.EQ.C.) EP*EE/1. 732051 B 108 

IF (EE.EQ.C.) EE=EP*l. 732051 B 109 

IF (GP.EQ.C.) GP=GC B 110 

IF (DWI.NE.G.) SH=CW1 B 111 

IF ^IPX.EG.0.A^C.RPM.^E.C. ) IPX=(F*120. »/RPM B 112 

PX=IPX B 113 

IF (RPM.EC.0..AN0.PX.NE.C. ) RPM= ( F«120. ) /PX B 114 

IF (F.EC.O.) F=PX*RPM/12C. B 115 

hW=HY-HC-HT B 116 

CC^IQQ B 117 

IF (ZZ.NE.3) GO TO 2 B 118 

ei=(H0+HT-FS)*{6.283ie5/CQ)+83 B 119 

e2=Bl+(6.283185*Hh/QQ) B 120 

8S=(B2+e3)/2. B 121 

2 CONTINUE B 122 

PI=(VA*10CC. )/«EE»SQRl (3.) ) B 123 

CK=l. B 124 

IF (PF.GE.C.95) CK=1.10 8 125 

IF (ZZ.E0.1.0R.ZZ.EQ.5) BO=BS B 126 

111=11 B 127 

CB=.25 B 128 

IF (DU.GE.8.> 08=0.5 B 129 

IF (BCOIL.EC.O. ) BCOIL=ALH B 130 

FE=3.1416*(PCCIL+CCOIL)/2. B 131 

CR=CI-2.*GC B 132 

IF (PE.EQ.O.) PE=<PX/3.1415927)»{ARSIN(PHW/DR)) 8 133 

IF (PHW.EC.O.) PHN=DR*SIN(3.1415927*PE/PX) B 134 

IF (BP.EO.C.) BP=PHW B 135 

IF (PL.EQ.C.) PL=PEL B 136 

IF (PHL.EC.O.) PHL=PL B 137 

PC=(0U-CI-2.C«HS )«C.5 B 138 

IF (OYC.EG.O.) CYC=OU B 139 

ZY=0.7*FS B 140 

DO 3 1=1.5 B 141 

3 IF (G(I).GT.9.) G(I)=G(I )/10C. B 142 

CN=CG/(PX*PN) B 143 

CS=YY/(PN*CN) 8 144 

B 145 

CHECK FOR ERROR CONDITIONS B 146 

B 147 

IF (CS.GT.l.C. CR.CS. lt. 0.5) WRITE 16,5) CS B 148 

IF IEP*EE.EC.O..OR.ABS(Et/EP-l.732C51).GT.O.Ol) WRITE (6,6) B 149 

IF (PX*F*RPM.EQ.0..OR.ABS(F-PX*RPM/12O.).GT.C.l) WRITE (6,7) B 150 

IF (HC.LT.ZY) WRITE (6,8) HC,HS B 151 

IF (DSH.GE.CR) WRITE (6,9) B 152 

IF (DCOIL.GT.DYC) WRITE (6,10) B 153 

IF (PC0IL.LT.Dl+2.«HS ) WRITE (6,11) B 154 

IF (TYPY.GT. l.ANC.TYE*TYR.LT.1.0E-10) WRITE (6,12) B 155 

IF (RT.LT. l.OE-10) GO TO 4 B 156 

IF ( ( ( (CCCIL-PC0IL)*6CCIL)/(RT»R0) ).LE.2.*PT) WRITE (6,13) B 157 

GC TO 14 B 158 

IF ((0CC1L-PCCIL)»BC0IL/RD«*2.LE.1.7146»PT) WRITE (6,13) B 159 

FORPAT (5X,27H CS (PER UNIT POLE PITCH) =,F7.3/10X,31H CS MUST BE B 160 

IBETWEEN 0.5 AND 1.0) B 161 

6 FORMAT (IH ,38H EITHER PHASE CR LINE VOLTAGE IS WRONG) B 162 


50 


SINDUC 



7 

8 

9 

10 

11 

12 

13 

C 

C 

C 

14 


15 


16 


17 


18 

19 

20 

21 

C 

c 

c 


22 

23 

24 

25 

26 

27 

28 
29 


FORMAT (IH ,44H FREQUENCY* RPM* CR NO. OF POLES IS IN ERROR) 

FCRPAT (IH /5X54HCEPTF BELCW SLOT IS LESS THAN 70 PERCENT OF SLOT 
1CEPTH/10X*4HDBS=F8.4/1CX,4H S0=F8.4) 

FORMAT (IH *46H SHAFT DIAMETER IS GREATER THAN ROTOR DIAMETER) 
FORMAT (IH *34H FIELD COIL O.D. EXCEEDS YOKE I.O.) 

FORMAT (IH ,29H FIELD COIL I.D. IS TOO SMALL) 

FORMAT (IH ,49H TYE AND TYR MUST BE READ IN FOR TYPE 2 OR 3 YOKE) 
FORMAT (IH #81H FIELD COIL DIMENSIONS ARE TOO SMALL FOR THE SPECIF 
IIEC NO. OF TURNS AND kIRE SIZE) 

DETERMINE ROTOR AND STATCR STACKING FACTORS 


M=1 

STFK^SF 
LT=LTS 
GO TO 17 
M==2 

STFK=RK 
LT=LTR 
GO TO 17 
M = 3 

STFK=RK1 

LT=LTR1 

IF (STFK.NE.O.) GO TO (19, 20, 21), M 
IF (LT.EQ.O.) GC TC 18 
STFK = l.C-( 12.5E-4/LT) 

GO ro (19, 20, 21), M 
STFK=1.0 

GC TO (19, 20, 21), M 

SF=STFK 

GC TO 15 

RK=STFK 

GO TO 16 

RK1=STFK 

CALCULATE POLE FACE LOSS FACTOR 


M=0 

IF (Dl.NE.O.) GO TC 29 
IF (LTRl.NE.C.) GC TO 22 
M=1 

IF (RKl.GT. 0.9999) GO TC 28 
LTR1=( 12.5E-4)/( 1.0-RKl) 

IF (LTR1-C.045) 23,23,24 

01=1.17 

GO TO 29 

IF (LTRl-0.094) 25,25,26 

01=1.75 

GC TO 29 

IF (LTRl-0.17) 27,27,28 

Dl=3.5 

GC TO 29 

Cl=7.0 

IF (M.EC.l) LTR1=0. 

IBN=8N+.l 

SS=SF*(CL-HV*BV) 


B 163 
B 164 
B 165 
B 166 
B 167 
B 168 
B 169 
B 170 
B 171 
8 172 
B 173 
B 174 
B 175 
B 176 
B 177 
B 178 
B 179 
B 180 
B 181 
B 182 
B 183 
B 184 
8 185 
B 186 
8 187 
B 188 
B 189 
8 190 
B 191 
8 192 
B 193 
6 194 
B 195 
B 196 
B 197 
B 198 
B 199 
B 200 
8 201 
B 202 
B 203 
B 204 
B 205 
B 206 
B 207 
B 208 
B 209 
B 210 
B 211 
B 212 
B 213 
B 214 
B 215 
B 216 
B 217 
B 218 
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SIGf'A==(54.E3/CI**2)*(PF/SS)*(VA/RPM) B 219 

VR=0.262*CR*RPM B 220 

TS=3.142*DI/CG B 221 

IF (2Z-4) 30,31f30 B 222 

30 1T=(.667*HS+DI)*3.142/CG B 223 

GC TO 32 B 224 

31 TT=(0I*2.0«HC+1.333*BS)*3.1416/GQ B 225 

C B 226 

C CALCULATE CARTER CCEFFICIENTS B 227 

C B 228 

32 IF (ZZ.GT.1.AND.ZZ.LT.5) GO TO 33 B 229 

CC=(5.0*GC+BS»*TS/n5.C*GC+BS)*TS-BS*BS) B 230 

GO TO 34 6 231 

33 CC=(4.44*GC+C.75»B0)*TS B 232 

CC=CC/(QC-BC*BO) B 233 

34 IF (IBN.EC.O) GC TO 35 B 234 

CC=(4.44*GC+0.75»W0)«TB B 235 

CCR=QC/(QC-hC«*2) B 236 

GC TO 36 B 237 

35 CCR=1. B 238 

36 TP=3.142*CI/PX B 239 

C B 240 

C PITCH FACTOR AND SKEW FACTOR CALCULATIONS B 241 

C B 242 

CF=SIN(YY*1.571/(PN*0M) B 243 

IF (SK) 37,37,38 B 244 

37 FS=1.0 B 245 

GC TO 39 B 246 

38 FS=(SK/TP)*1.5707 B 247 

FS=(l./FS)*(SIN(FS))*CCCS(FS*(l.+BCOIL/CLm B 248 

C B 249 

C CHECK IF WINCING HAS INTEGRAL NO. OF SLOTS PER POLE PER PHASE B 250 

C B 251 

39 0*1.0 B 252 

IF (PBA.GT.61.0) C=2.C 8 253 

IZY*IPX»IPN B 254 

ICH-=0 6 255 

40 IDH=IDH+IZY B 256 

IF (IQC-ICH) 42,41,40 B 257 

C B 258 

C CALCULATE CISTRIBUTION FACTOR FOR INTEGRAL SLOT WINDING 6 259 

C B 260 

41 CF*SIN(1.571»D/PN)/(CN«D»SIN(1.571/(PN*QNJ n B261 

GC TO 46 B 262 

C B 263 

C CALCULATE DISTRIBUTION FACTOR FOR FRACTIONAL SLOT WINDING 8 264 

C B 265 

42 IICC=ICC 6 266 

1*2 B 267 

43 IF ( (IZY/n«I.EQ.IZY.ANO.( I1C0/I)*I.EQ.IIQQ) GO TO 44 B 268 

IF (I.GT.IZY) GC TC 45 B 269 

1*1*1 B 270 

GC TO 43 B 271 

44 IZY=IZY/I B 272 

IICC=IICQ/I B 273 

GO TO 43 B 274 
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45 FKC=°IICC B 275 

CF=SIN(1.571*C/PN»/(FNC«C*SIN(1.571/(FNQ*PN))) B 276 

46 EC=CO*SC*CF*FS/C B 277 

C B 278 

C CCPFUTE ABPATURE CCNOLCTCR AREA B 279 

C B 280 

IF (OMl) 47*47*48 B 281 

47 AC=C.765*0k«0h»SNl B 282 

6C TO 60 B 283 

48 ZY=0.0 B 284 

OT»AMlNl(Ck*DMl) B 285 

C6=AMAXl(CW*Ckl) B 286 

49 IF (0T-.05I 52*52*50 B 287 

50 JA=0 B 288 

51 JA=JA+1 B 289 

IF (OT-DAUA)) 53.53,51 B 290 

52 C=0 B 291 

IF (ZY» 59,59,72 B 292 

53 IF (DG-0.188) 54,54,55 B 293 

54 CY-CX(JA-l) B 294 

CZ*CXIJA) B 295 

GO TO 58 B 296 

55 IF (DG-0.75) 56,56,57 B 297 

56 CY = CYUA-1> B 298 

CZ=CYIJA) B 299 

GO TO 58 B 300 

57 CY=CZCJA-1) B 301 

CZ-CZIJA) B 302 

58 C=CY+(CZ-CY)«(DT-CA(JA-in/(CA(JA)-DA(JA-l)) B 303 

IF (ZY) 59*59*72 B 304 

59 AC=(DT*CG-D)*SN1 B 305 

C B 306 

C CALCULATE END EXTENSICN LENGTH B 307 

C B 308 

60 IF (EL) 61,61,69 B 309 

61 IF (RF) 62*62,68 B 310 

62 IF (PX-2.0) 63,63,64 8 311 

63 L=1.3 B 312 

GO TC 67 B 313 

64 IF (PX-4.0) 65,65,66 3 314 

65 L=1.5 6 315 

GO TO 67 6 316 

66 U*1.7 B 317 

67 EL=3.1416»U«YY*CDI*HS )/(3C+C.5 B 318 

GO TO 69 8 319 

68 EL=2.0*CE+(3.1416*(0.5«HX+DBn + (YY*TS*TS/(SQRT(TS*TS-BS*BS) n 6 320 

69 HN=2.*CL+EL+BCCIL B 321 

C 8 322 

C CALCULATE STATOR RESISTANCE B 323 

C B 324 

A*PI*SC*CF/(C*TS) B 325 

RY=SC*CO*HP/(PN*AC«C*C ) B 326 

RGl=(l.E-6)*RS*(l.C+ALPHAS«(TST-20. n*RY B 327 

S*PI/(C*AC) 8 328 

C 8 329 

C COMPUTE FIELD CONDUCTOR AREA B 330 
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c 


B 

331 


IF (RT) 70,70,71 

B 

332 

70 

AS=.7854*RC»RD 

B 

333 


GC TC 73 

B 

334 

71 

ZY^rl.O 

B 

335 


DT=/^MINl(RT,RD) 

B 

336 


DG = /iMAXnRT,RO) 

B 

337 


GC TO 49 

B 

338 

72 

A$:sCT*CG-C 

B 

339 

C 


B 

340 

C 

CGf'PUTE FIELD RESISTANCE 

B 

341 

C 


B 

342 

73 

2G=PT*FE/AS 

B 

343 


FKl=( l.E^6)*RR*( 1.0+ALPHAR»(TF-2C. ) )*ZG 

B 

344 

C 


B 

345 

C 

NC LOAO MAGNETIC CALCLLATICNS 

B 

346 

C 


B 

347 


G A* 3. 142*01* (CL-HV*BV ) 

B 

348 


GE*CC*GC*CCR 

B 

349 


AG=6.38*0I/(PX*GE) 

B 

350 


IF (Cl) 75,74,75 

B 

351 

74 

C1=(.649*ALCG(PE)+1.359)*{ ( GC/GP ) * *0 . 352 ) 

B 

352 

75 

CU*C.70 7*EE»C1*CF/(EP*PN ) 

B 

353 


TG=6OOCC0C.0«EE/ (CV^*EC*RPM) 

6 

354 


BG=TG/GA 

B 

355 


IF (CP) 76,76,77 

B 

356 

76 

CP=(GC/GP)**.41*PE*( ALCG(GC/TP)*.0378+1.191) 

B 

357 

77 

FQ^TG*CP/PX 

6 

358 

C 


B 

359 

C 

DETERMINE DEMAGNETIZING AMPERE TURNS (FULL LCAO) 

B 

360 

C 


8 

361 


IF (CM) 78,78,79 

6 

362 

78 

AA==SIN(3*142*PE) 

B 

363 


AB*SIN( 1*571*PE)*4.0 

B 

364 


CM=(3.142*P£+AA)/AB 

B 

365 

79 

CONTINUE 

B 

366 


FGML=.45*EC*PI*CM*CF/PX 

B 

367 

C 


6 

368 

C 

PERMEANCE CALCULATIONS 

6 

369 

C 


B 

370 


IF (CQ) 80,80,61 

B 

371 

80 

AB=*3.1416*PE 

B 

372 


CG=(4.*PE41. )/5.-SIN( Ae)/3,1416 

6 

373 

81 

XR==.0707*A*CF/(C1*BG) 

B 

374 


FACTCR-YY/ (PN*CN) 

B 

375 


IF (PBA.LT.61.) GC TO 82 

6 

376 


FF==.05*(24.*FACTCR-“1, ) 

B 

377 


IF ( FACTOR. GE. 0*66 7) FF=.75 

B 

378 


IF (ZZ.EQ.5) FF=1. 

8 

379 


GC TC 83 

6 

380 

82 

FF«.25* (6.*FACTCR-1. ) 

B 

381 


IF (FACTOR. GE. 0.667) FF==.25* (3.*FACT0R+1. ) 

B 

382 


IF (ZZ.EQ.5) FF^^l. 

B 

383 

83 

CX=FF/(CF*CF*CF*DF) 

B 

384 


Z=CX*2C.0/(PN*CN) 

B 

385 


BT=:3.142*DI/CQ-B0 

B 

386 
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ZA=eT*BT/( 16.0*TS*GC ) B 387 

ZB=0,35»BT/TS B 388 

ZC=hC/0C B 389 

ZC=hX*.333/BS B 390 

ZE=hY/8S B 391 

IF (ZZ-2) 84,85,06 B 392 

84 PC=Z*(ZE+ZC+ZA+ZB) B 393 

GC TO 90 B 394 

85 PC = Z*(ZC+(2.0*HT/(BO+eS) )+(HW/BS)+ZD+ZA + ZB) B 395 

GO 10 90 B 396 

86 IF (ZZ-4) 87,88,89 6 397 

87 PC=Z*(ZC+(2.0*HT/(EO + E1) )+ ( 2 . 0*hV</ ( B 1+B2 ) ) + ( HX/ ( 3. *62 ) ) + ZA + ZB) 6 398 

GO TO 90 B 399 

88 PC=Z*(ZC+0.62) B 400 

GC TO 90 8 401 

89 PC=Z*(ZE+ZC+(0.5*GC/TS)+(0.25*TS/GC)+C.6) B 402 

90 EK=EL/( 10.0**(0.103*YY*TS+0*402) ) B 403 

IF (DI-8.C) 91,91,92 B 404 

91 EK=SQRT(EK) B 405 

92 ZF=.612*ALCG{10.0*CS) B 406 

EW=6.28*EK*ZF* (TP**( 0.62- ( .228* ALOG(ZF) ) ) ) / ( CL*DF*DF ) B 407 

PM=3.19*3.1416*CR*CL*(2.C-PE)/(PX*(HP1+GC) ) B 408 

P5=l.675*{CCCIL-PC0IL )*{ CCCIL + PCOID/BCOIL B 409 

P6 = 2.5*(PCCIL-DI )*(PCC IL+DI )/BCOIL B 410 

P6 = P6+1.67* (DI-DSH)* < C I + CSH ) / BCD I L B 411 

P7 = 2.5*(DI+DISH1)*{0U-CI )/(Ot-DlSHl) B 412 

RL= (P5 + P6 + P7/2. + PK*PX/4. ) 6 413 

STATET=CQ*SC*DF*CF/( 2.«PN*C) B 414 

C B 415 

C STATOR WINCING LEAKAGE AND ARMATURE REACTION REACTANCES B 416 

C B 417 

XL=XR*(2.*PC+EW) B 418 

XD=XR*AG*Cl*CM B 419 

XC=XR*CC*AG B 420 

C B 421 

C FIELD LEAKAGE REACTANCE, SELF INDUCTANCE AND TIME CONSTANT B 422 

C B 423 

XF=3.0E-06*3.1416*F* { S T A TET • * 2 ) *RL*P I /EP B 424 

SI=PT*PT* (PX*3.1416*CP*AG*CL/8.+RL )*1.E-C8 B 425 

TC=SI/FK1 B 426 

C B 427 

C SYNCHRONOUS AND TRANSIENT REACTANCES CALCULATIONS B 428 

C B 429 

XA=XL+XD B 430 

XB=XL+XC 6 431 

XU=XL+(XF*XD)/ (XF+XD) B 432 

C B 433 

C COMPUTE FRICTION AND WINCAGE B 434 

C B 435 

IF (WF-1.0) 94,93,94 B 436 

93 WF=CR**2.5*(RPN**1.5 )*PL«O.OC000252 B 437 

C B 438 

C WEIGHT CALCULATIONS B 439 

C B 440 

94 IF (ZZ-3) 95,96,95 B 441 

95 WI=( (DU+DI)*(DU-DI)*3.1416)/4. B 442 
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IF (ZZ.NE.4) hI=WI-QQ*(BS*HS-( (H0+0.5*HT)*{BS-B0) » ) B 443 

IF (ZZ.eO.4) hI*m-QC«(BS*BS«3. 1416/4. + HO*BO) 6 444 

6C TO 97 B 445 

96 Wl=(DU-hC)«3.1416*HC B 446 

WI=ViI*HS*nDI+2.*HS)*3.1416-CQ«B3) B 447 

hI*hI+CC*( (HC+0.5»HT)«(BS-BOn B 448 

97 hI*hI*0.566*SS B 449 

C B 450 

RC»0.321*PT*FE«AS B 451 

C B 452 

WC*.321*SC*CC»AC*hF B 453 

C B 454 

IF (TYPY.EC.l) 6C TO 98 B 455 

Vi YOKE*. 283*1 ( 3 . 14 16* ( C YC +TYE ) *T YE* ( BCOIL+2 .*TYR ) )♦( 3.1416*1 (DU+TY+ B 456 

lCYC>/2.)*2.*TYR*(CYC-(CU + TYn/2.)) B 457 

IF (TYPY.EC.2) WY0KE=hYGKE+3. 1416*0. 283*(DU+TY)*TY*(2.*CU B 458 

IF (TYPY.EC.3) WYCKE=KYCKE+3. 1416*0. 283*2. *CL*(0. 333*1 10.5*DU+TY)* B 459 

1*2 + 1 0.5» lot ♦TYJ )**2*(0.5*10U+TY> )*10.5*DL+TY) )-0 .25*DU*L)U ) B 460 

GO TO 99 B 461 

98 WYOKE*. 283*3. 1416*1DU+TY)*TY*12.*CL+BC01L) B 462 

C B 463 

99 IF IWRCTOR.NE.O. ) GO TO 100 B 464 

ViSHAFT*. 283*3. 1416*1DSF**2-D1SH**2)/4.*1ALH+2.*PL) B 465 

THETA=2.*3.1416*PE/PX B 466 

ATIP*DR**2*lTHETA-SINlTHETA))/8. B 467 

ABCCY*CR*SIMTHETA/2. )*lCR*CCSlTHETA/2. )/2.-CSH/2. ) B 468 

EETA*ARSIM 10R*SINlTHETA/2. ) /2. )/lDSH/2. n*2. B 469 

ABASE=0SH**2*lSINlBETA/2. )-SlNlBETA)/4. -BETA/4. )/2. B 470 

WPOLE=.283*PL*1ATIP+AECOY+ABASE) B 471 

V»ROTOR*ViSFAFT + PX*WPOLE B 472 

100 V«T0TAL*V»C + Wl + RC+HY0KE4hRCTCR B 473 

RETURN B 474 

END B 475- 
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SUBROUTINE MAGNET C I 

COMMON A, AA,AB,ACf ACRtAGtAI,ALH,ALPHAE,ALPHARf ALPHAStALY.ALYCtALYR C 2 

ItAPtASi ASH. ATH. AY, AYC* AYR. 8.61. B2tB3, BOLL .BC0 IL.bg tBKt BN, BO. BP. BPL C 3 

2,BS.B$HL.BTLL.BV.BYCLL,C.C1.CC.CCR.CE.CE,CK.CL,CM,CP,CO,CM,O.DI,DC C 4 

3CIL,OD.CF.OI.OISH.DISt-l.CR.OSH.OU.OW,OWl,OYC,EC.EOO.EE,EL,EP,EM.F, C 5 

4FCL,FE,FFLL,F6L,FGML, FF.FKl.FPL.FU.FS.FSHL.FSHLP.FTL.FYCL.FYL.FYRL C 6 

5,G,GA,GC,GE,GP,GXX,H,HC,H0,HM,H0,HP,HP1,HS,HT,HV,HW,HX,HY,I8N,IPN, C 7 

6IPX,IQQ,IZZ, JA,KSAT,LTR,LTRl,LrS,P5,P6,P7,PBA,PC,PC0IL,PE,PF,PHU,P C 8 


7HW,PI ,Pt,PP,Pf^L,PN,PT.PX»QN,CQ,RC 

,RD,RE,RFfRGl,RK,RKl,RPP,RR,RS,RT 

C 

9 

8,RY,S,$B,SC,SC,SF,SH,SI,SIGMA,SK, 

SN,SN1,SS,STATET,T1,T11,T2,T22,T3 

C 

10 

9*T33,TE,TC,TF,TG,TS,TST,TT,TY,TYE 

,TYPY,TYR,VA»VR#WC.WF,WI ,WL,WO, WR 

C 

11 

$CTCR,WTGTAL,KYCKE,XA,>B,XDfXF,XL, 

XQ,XR,XU,YY,Z»ZG,ZZ,ZZZ 

C 

12 

INTEGER TYPY 


C 

13 

CIKENSICN /1I(90), G(5) 


C 

lA 

BPL = 0 


C 

15 

BTLL=0 


C 

16 

BSHL=0 


C 

17 

BCLL=0 


c 

18 

BYCLL=0 


c 

19 

FFLL=0 


c 

20 

FTL = 0 


c 

21 

FPL = 0 


c 

22 

PPL=GXX*FC 


c 

23 

W=FH*ECC 


c 

2A 

FGL^W 


c 

25 

PHLA=0 


c 

26 



c 

27 

CALCULATE AIRGAP Af^PERE-TURNS ANC 

PML (LEAKAGE FLUX BTWN POLES) 

c 

28 



c 

29 

P/«L = PM* (FGI^L + FTL + FPL + FCL )*.OCl 


c 

30 

P^«LAG^P^^L^lPE/(2.0-PE ) 


c 

31 

FGL=W+PRLAG*ZZZ 


c 

32 

PRL=PM* (FGNL+FTL+FPL+FGL)*.CC1 


c 

33 

PMLAG^PPL«PE/(2.0-PE) 


c 

3A 



c 

35 

FLLX DENSITY ANC A^P E R E-- TURNS FOR 

POLE 

c 

36 



c 

37 

BPL=(PPL+PRLAG)/AP 


c 

38 

NA = 31 


c 

39 

K = 1 


c 

AO 

X = BPL 


c 

A1 

GO TO 19 


c 

A2 

FPL=AT*HP 


c 

A3 



c 

AA 

FLUX DENSITY AND AHPERE-TURNS FOR 

TEETH 

c 

A5 



c 

A6 

BTLL=(PPL+PRLAG)/ATH 


c 

A7 

X=BTLL 


c 

A8 

NA=l 


c 

A9 

K = 2 


c 

50 
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GO TO IS 
3 FTL=AT*HS 

C 

C CHECK IF Pf^L HAS CONVERGED 

C 

IF (ABS((PHL-PRLA)/PHL).LE.l*CE-04) GO TO 4 
PRLA=PRL 
GO TO 1 
C 

C FLUX DENSITY AND AHPERE--TURNS FOR SHAFT (UNDER FIELD COIL) 

C 

A Z=FTL+FGL+FPL 

PH7L=P7*Z*.001 

PSHL=(PPL+PRLAG)*PX/2.C+PML*PX/2.0+PH7L 

BSHL=^PSHL/ASH 

X^BSHL 

NA=31 

K==3 

GO TO IS 

5 FSHL=AT*ALH 
C 

C FLUX DENSITY AND Af^P E P E-TURNS FOR SHAFT (UNDER POLES) 

C 

PDIFF=PSHL-PH7L 
X= ( .250*PC IFF + PH7L )/ASH 
NA = 31 
K = A 

GO TO IS 

6 FSHLP=AT*PL/2«0 
X=(.625*PCIFF+PH7L)/ASh 
NA = 3l 

K=5 

GO TO IS 

7 FSHLP=FSHLP+AT»PL/A.O 
X= ( .875*PCIFF + PH7L )/ASH 
NA = 31 

K = 6 

GO TO IS 

8 FSHLP=FSHLP + AT^^PL/A.0 
C 

C FLUX DENSITY AND ARPEPE-TURNS FOR CORE 

C 

Z=2.*Z+FSHL+FSHLP*2. 

PH6L=P6*Z*.C0i 

eCLL=(PPL + PPLAG+( PF7L4PH6L )/PX )/ACR 
X=BCLL 
NA=1 
K = 7 

GC TO IS 
S FCL=AT*HC 

C 

C FLUX DENSITY AND APPEPE-TURNS FOR YOKE 

C 

Z=^Z + 2.*FCL 

PH5L=P5*Z*.C01 

IF (TYPY-1) 11,10,11 


C 51 
C 52 
C 53 
C 5A 
C 55 
C 56 
C 57 
C 58 
C 5S 
C 60 
C 61 
C 62 
C 63 
C 6A 
C 65 
C 66 
C 67 
C 68 
C 6S 
C 70 
C 71 
C 72 
C 73 
C 7A 
C 75 
C 76 
C 77 
C 78 
C 7S 
C 80 
C 81 
C 82 
C 83 
C 8A 
C 85 
C 86 
C 87 
C 88 
C 89 
C 90 
C 91 
C 92 
C 93 
C 9A 
C 95 
C 96 
C 97 
C 98 
C 99 
C 100 
C 101 
C 102 
C 103 
C lOA 
C 105 
C 106 
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I 



10 

PY=PSHL+PH6L+PH5L 

C 

107 


GC TO 12 

C 

108 

11 

PY=PSHL+PH6L 

C 

109 

12 

X-PY/AY 

C 

110 


NA = 61 

C 

ill 


K = 8 

C 

112 


GC TO 19 

C 

113 

13 

FYL=AT»ALY 

C 

114 


IF (TYPY-1) 14fl5,14 

C 

115 

14 

PY-PY+PF5L 

C 

116 


X*PY/AYC 

C 

117 


BYCLL=X 

C 

118 


NA = 61 

C 

119 


K*9 

C 

120 


GO TO 19 

c 

121 

15 

FYCL==0 

c 

122 


FYRL*0 

c 

123 


BYCLL=X 

c 

124 


GC TO 18 

c 

125 

16 

FYCL=AT*ALYC 

c 

126 


X=PY/AYR 

c 

127 


NA = 61 

c 

128 


K=1C 

c 

129 


GC TO 19 

c 

130 

17 

FYRL-AT*ALYR 

c 

131 

C 


c 

132 

C 

TOTAL APPERE-TURNS 

c 

133 

C 


c 

134 

18 

FFLL*2.*(FGL+FTL+FCL+FPL+FSHLP )+FSHL+FYL+FYCL+FYRL 

c 

135 


RETLRN 

c 

136 

C 


c 

137 

C 

INTERPOLATION PROCEDURE FOR PATERIAL CURVES 

c 

138 

C 


c 

139 

19 

IF CAI(NA)-X) 24,2C,2C 

c 

140 

20 

NA=NA+3 

c 

141 

21 

IF (AI(NA)-X) 22f23,23 

c 

142 

22 

NA-NA+2 

c 

143 


GC TO 21 

c 

144 

23 

AA^AKNA) 

c 

145 


EB1=AI (NA-2) 

c 

146 


CC=AI(NA+1) 

c 

147 


C*AI(NA-1) 

c 

148 


XX=(AA-BBl)/(.4343*( ALCG(DC)-AL0G(D+.0001) )) 

c 

149 


Y=AA-XX*.4343*ALCG(0C ) 

c 

150 


AT=EXP(2*3C6«(X-Y)/XX) 

c 

151 


GO TO (2,3»5,6»7,8»9, 13, 16,17) ,K 

c 

152 

24 

KSAT=0 

c 

153 


RETURN 

c 

154 


END 

c 

155- 
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SUBROUTINE OUTPUT 

COFPON A, AA,A6f AC,ACk,/16,AI ,ALH,ALPHAE,ALPHAR,ALPHAS,ALY,ALYC»ALYR 0 2 

ItAP.ASt ASH»ATHt AY, AYCt AYRtBtBl ,B2>B3,BCLL,BC0ILfBG«BK,BN,B0,BP,BPL 0 3 

2,ES,BSHL,BTLL,BV,eYCLL,C,Cl,CC,CCR,CE,CF,CK,CL,CM,CP,CQ,CW,D,DI,DC D 4 

3CIL,DO,CF,CI,DISH,ClSH,CR,DSH,DU,DW,DWl,CYC,EC,EDD,EE,EL,eP,EW,F, D 5 

4FCL,FE,FFLL,FGL,FGPL,FH,FKl,FPL,FQ,FS,FSHL,FSHLP,FTL,FYCL.FYL,FYRL D 6 

5,G,GA,GC,GE,GP,GXX,H,HC,FD,HP,H0,HP,HP1,HS,HT,HV,HW,HX,HY, IBN, IPN, D 7 

6IPX,IQQ,IZZ, JA,KSAT,LTR,LTRl,LTS,P5,P6.P7,PBA,PC,PC0IL,PE,PF,PHL ,P 0 8 

7HW,PI,PL,PP,PNL,PN,PT,PX,QN,CQ,RC,RD,RE,RF,RGl,RK,RKl,RPP,RR,RS,RT D 9 

8, RY,S,SB,SC,S0,SF,SH,SI,SIGHA,SK,SN,SNl,SS,STArET,Tl,Tll,T2,T22,T3 D 10 

9, T33,TB,TC,TF,TG,TS,TST,TT,TY,TYE,TYPY,TYR,VA, VR,WC,WF,WI,WL,WO,WR 0 11 


$CTCR,WTOTAL,ViYCKE,XA,xe,XO,XF,XL,XQ,XR.XU,YY,Z,ZG,ZZ,ZZZ 0 12 

C 0 13 

DIMENSION STAR(5), DASH(3), AI(90). G(5) D 14 

C 0 15 

INTEGER TYPY 0 16 

C 0 17 

real LTS.LTR.LTRl 0 18 

C 0 19 

DATA STAR ( I ) /30H«««** ••••*•***•••••**••*******/, OASH ( 1 )/18H 0 20 

1 / 0 21 

0 22 
0 23 

WRITE 16,1) VA,EE,EP,PI,FF,IPN,F,IPX,RPH 0 24 

FORMAT (1HL,18H ALTERNATOR RATING//10X, 15H ALTERNATOR KVA, F16. 1/10 0 25 

1X,18H LINE-LINE VOLT AGE , F12.0/10X, 19H LINE-NEUT. VOLTAGE , F 1 1 .O/lOX 0 26 

2,14H PHASE CURRENT, Fie. 2/10X,13H POWER F ACTOR , F19. 2/ lOX, 7H PHASES, 0 27 

3I22/10X,1CH FRECUENCY,F2C.0/10X,6H POLES, I23/10X,4H RPM,F27.1) 0 28 

IF (IZZ-2) 3,5,2 0 29 

IF (IZZ-4) 7,9,11 0 30 

WRITE (6,4) BS,HX,HY,FS, IQC,TS,TT 0 31 

FORMAT (1HL,13H STATOR SLOTS//5X10H TYPE-0PEN/54X, 9H »,12X6 0 32 

IH* /62XlH*,12XlH*/55X2HHY,5Xl(i*, 12XIH*/10X3H BS , F 26 . 3 , 1X6H INCH 0 33 

2ES, 16X, 1H*,12X1H»/10X3H hX , F26 . 3 , 15X , 9H * , 2X8H»*«**»» * , 2X IH D 34 


3«/lCX3H HY,F26.3,23XIM,2XIH*,6X1H*,2XIH»/10X3H HS , F26 . 3 , 23X IH* , 2X 0 35 

41H* ,6X1H»,2X1H*/62X, IF •,2X8H« »*••***, 2X1H*2X 2HHS/55X2HHX,5X, IH*, 12 0 36 

5X1H*/10X13H NC, OF S LCTS I 16 , 23X , IF* , 2X8H** ** •*** , 2X IH* /62X1H* , 2X IH 0 37 

6*,6X1H* ,2X1H*/10X11H SLOT P I TCH , F 1 8. 3 , 1X6H INCHES , 16X IH* , 2X IH* , 6X IH 0 38 

7*,2X1H*/54X9H • ,2X8H********,2XIH*/ 10X1 IH SLOT PITCH, 41XIH* 0 39 

8,12X1H*/10X15H AT 1/3 C I S T . , F 14. 3, IX6H INCHES , 16X 19H************* 0 40 

9* /62X1H1,12X1H1/62X14H1 BS 1/62X1H1 , 12X1H I ) 0 41 

GO TO 13 0 42 

5 WRITE (6,6) BC,BS,HO,FX,FT,HW,HS,ICQ,TS,TT 0 43 

6 FORMAT (1HL,13H STATOR SLOTS//5X22H TYPE-PART I ALL Y CL0SE0/67X4H-B0 0 44 

1-/57X10H *,4X1CH* /58X2HH0 , 6X IH* , 4X IH* /5 7X lOH 0 45 

2 *,4X1H*/10X3H BO,F26.3,1X6HINCHES,19XIH*,6X1H*/10X3H BS,F26.3, 0 46 

319X2HHT,4X1H*,8X1H*/1CX3F HO , F26 . 3 , 24X1H*, lOX IH*/ 10X3H HX,F26.3,18 0 47 

4X6H *,12X1H*/10X3F HT,F26.3,23XIH*, 12X1H*/10X3H H W , F26 . 3 , 19X2H 0 48 

5FW,2X1H*, 12X1H*/10X3H HS ,F26 . 3 , 18X6H *2X8H**** **** ,2X IH* , 2X2HH 0 49 

6S/62X1H*,2XIH*,6XIH*,2XIF*/1CX13F NO. OF SLOTS 1 16 , 23X IH* , 2 X IH* ,6X 1 0 50 
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7H*,2XlH*/62XlH*,2X8H*«««<HMi*t2XlH*/l0XllH SLOT PITCH,F18.3, 1X6HINC 
8HES,12X2HHX,2XlH*»12XlF*/62XlH»,2X8H********f2XlH*/10XlXH SLOT PIT 
9Ch,AlXlh*,2XlH*,6XlH*,2XlH*/10X15h AT 1/3 Cl ST, ,F14.3, 1X6HINCHES 

$t 16XIH*2X1H*6X1H*2X1H*/57X6H 2X8H*** **«**♦ 2X IH*/62X IH* • 12X IH 

$*/62X19H* **«•*<* «««♦*« /62XIH1, I2X1H 1/62X1 4H I BS 1/62X 

$1H1.12X1H1) 

GC TO 13 

7 V^RITE (6,8) BO,Bl,B2,e3.BS,HO,HX,HT,HW,HS,IQQ,TS,TT 

8 FORPAT (1HL,13H STATOR SLOTS//5X25HTYPE-CONSTANT TOOTH WI0TH/61X1H 

11, 14X1H1/61X16H1 B1 1/10X3H BO , F26 . 3 , 1X6H INCHE S , 15X IHl , 1 

24X1H1/1CX3H ei,F26,3,22XIHl,5X^H-BO-,5XlHl/lCX3H B2 , F26, 3 , 1 IXl 7H — 

3 1 #,4X17H» 1 /10X3H B3, F26, 3, 22X1H1 ,4X1H* 

4,4X1H*,4X1H1,8X2FH0/1CX15H BS * ( B2^^B3 ) /2, F 14 . 3, 22X1H1 ,4XIH* ,4X17H 

5, 1 /10X3H H0,F26,3,22X1HI,2X1H*,8X1H*,2X1H1,8X2HHT/1 

60X3F HX,F26.3,22X1F*, 14X,12H* /10X3H HT ,F^6 . 3 , 12X2HHS, 7X 

71H*,16X1H*,7X2HHW/10X3H HW , F26, 3, 20X IH* , 3X 12H***»* ***»*♦*, 3X10H* — 

8 /10X3H HS,F26. 3,19X2H«1,3X1H*,10X1H*,3X2H1*/57XIH*,1X1H1,3X 

91H*,10X1H*,3X1H1,1X1H«,4X2HHX/10X13H NO, OF SLOTS, I 16 , 17X IH* ,2X1H1 

$, 3X12H« ********** *,3XIH I, 2X1H*,6H /55X IH* , 3X IH 1 , 18X IH 1 , 3X1H#/ 

SlOXllH SLOT PITCH,F18*3, 1X6H INCHES ,4X34H ********************** 

$«*»«**«*/ 54X IH I ,4X1H1,18X1H1,4X1H1/ 10X1 IH SLOT P I TCH , 33X IH I , 4X20H 1 

$ B2 1,4X1H1/10X15H AT 1/3 01 ST • , F 14 . 3 , IX6H INCHES , 8 

$X1HI,4XIH1,18XIH1,4X1F1/54X3CH1 B3 1/54X1H 

$1,28X1H1) 

GC TO 13 

9 WRITE (6,10) BC,hC,BS,FS,IGQ,TS,TT 

10 FORPAT (1HL,13H STATOR SL0TS//5X, I IH TYPE-ROUND// lOX, 13H SLOT OPEN 
1ING,F16.3,1X6HINCHES/ICX,19H SLOT OPENING DEPTH, F 10 • 3/ lOX , 14H SLOT 
2 DIAMETER, F15.3/10X11F SLOT CEP TH , F 18 • 3// lOX , 1 3H NO. OF SLOTS, 116/ 
3/lOX,llH SLOT PITCH,F18.3, 1X6HINCHES//10X, IIH SLOT P I TCH/ lOX , 15H 

4 AT 1/3 DIST.,F14.3, 1X6HINCHES) 

GO TO 13 

11 WRITE (6,12) BS,HX,HY,FS.ICQ,TS,TT 

12 FORMAT (1HL,13H STATOR SLOTS//5X25H TYPE-OPEN (1 COND. /SLOT ) /57X ,6 

Ih *12X6h* /62X, IF*, 12XIH*/58X5HHY * , 1 2X 1H*/62X IH* , 12X IH*/ 

210X,3H BS,F26.3, 1X6H INCHES, 1 IX, 6H *, 2X8H*« ****** ,2X1H*/ 10X,3H 

3HX,F26.3,23X, 1H*,2XIH*,6X1H*,2X1H*/10X,3H HY , F 26 . 3 , 2 3X , IH* , 2X IH* , 6 
4Xlh»,2XlH*/10X,3H FS , F 26 .3 , 2 3X , IH* , 2 X IH* ,6X1H* ,2X IH* , 2X2HH S/58X2HH 
5X,2X1H*,2X1H»,6X1H*,2X1H*/10X,13H NO. OF SLO TS , 1 16 , 2 3X IH* , 2X IH* , 6X 
61H* ,2X1H*/62X1H*,2X1H* ,6Xlh*,2XlH*/10X, IIH SLOT PI TCH, FI 8. 3, IX6H1N 

7CHES,16XIH*,2XIH*,6X1F«,2X1H*/57X6H *, 2X8H*** »****, 2X IH*/ 10X1 1 

8H SLOT PITCH, 41X1F*, 12X1F*/1CX15H AT 1/3 D I ST . , F 14 . 3 , 1X6H I NCHES , 

9 16X19H* ************* /62XIH1,12X1H1/62X14H1 BS 1/62X1H 

$1,12X1H1) 

13 CONTINUE 

WRITE (6,14) GC,GP,GE 

14 FORMAT (1HL,8H AIR GAP//10X,16H MINIMUM AIR GAP , F 17. 3, IX6H INCHES/ 1 
ICX,16H MAXIMUM AIR GAP , F 17. 3/ lOX, 18H EFFECTIVE AIR GAP,F15.3//) 

IF (IBN.EC.O) GC TO 16 
WRITE (6,15) CC,CCR 

15 FORMAT (IH , lOX , 18HC AR TER COEFF IC lENT/ 17X , 6HST ATOR , F20. 3/ 18X5HR0T0 
1R,F20.3) 

GC TO 18 

16 WRITE (6,17) CC 

17 FORMAT (IH , lOX , 18HCARTER COEFF IC I ENT , F14.3 ) 

18 IF (RF.LT..5) WRITE (6,19) 
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IF (RF.GE..5) hRITE (6,20 D 107 

19 FCRNAT (1H1,45H ARMATURE WINCING ( Y-CONNECTED , RANDOM WOUND)//) 0 108 

20 FORMAT (1H1.A3H ARMATURE WINCING ( Y-CONNECTEO , FORM WOUND)//) 0 109 

IF (DWl.EC.O.) WRITE (6,21) CW D 110 

21 FORMAT (IH ,9X,16H STRAND D I AME I ER , F32 . 4 , 1 X6HI NCHES ) 0 111 

IF (DWl.GT.O.) WRITE (6,22) CW,DW1,SH 0 112 

22 FORMAT (IH ,9X,18H STRAND DIMENSIONS, F30. 4, 2H X, 1XF6.A, 1X6HINCHES/ D 113 

110X,35H UNINSULATED STRAND HEIGHT (RADIAL) ,F13. 4) D 114 

WRITE (6,23) SD,SN,SN1,AC,S,CE,HM,EL,SK,RS,TST,RG1,STATET,YY,QN D 115 


23 FORMAT (IH , 10X36HCI STANCE BTWN CL OF STRANDS (RAO I AL ) ,F 1 1.4//10X, D 116 
133H STRANDS/CONOUCTOR IN RADIAL D IR. , F 1 1 .O/lOX , 24H TOTAL STRANOS/C D 117 
2CN0LCTOR,F20.O/10X,15H CCNCUCTCR AREA,F33i4, 1X6HSQ-IN./10X,29H CUR D 118 
3RENT DENSITY AT FULL LC AC. F 17.2 . 3X lOHAMP/SQ- IN .//lOX ,27H COIL EXTE D 119 
4NSICN BEYOND CORE , F20 . 3, 2X6H INCHES/ 10X,24H MEAN LENGTH OF 1/2 TURN 0 120 
5,F23.3/10X,16H END TURN LENGTH, F31 .3/ lOX , 30H STATOR SLOT SKEW (PER 0 121 
6 STATOR), F17.3//10X,25H RESISTIVITY AT 20 DEG. C,F23.4, 1X16HMICR0 D 122 
7CHM INCHES/10X,21H STATOR RESISTANCE AT,F6.0,7H DEG. C ,F 14.4 , 1X4H0 D 123 
8HMS//10X,30H NO. OF EFFECTIVE SERIES TURNS , F 16 .2/lOX , 14H SLOTS SPA D 124 


gNNEC,F30.0/10X,25h SLOTS PER POLE PER PHASE,F21.2) D 125 

WRITE (6,24) SC,C,FBA,FS,DF,CF D 126 

24 FORMAT (IH ,9X,16H CONDUCTOR S/ SLC T , F2 8 . 0/ lOX , 2 5H NO. OF PARALLEL C 0 127 

1IRCLITS,F19.C/10X,17H PHASE BELT ANGLE , F27 . 0 , 5 X7HDE6REES/ / lOX , 12H D 128 
2SKEW FACTOR, F35.3/10X,20H DISTRIBUTION FACTOR, F27. 3/lOX, 13H PITCH 0 129 
3FACT0R,F34.3) D 130 

IF (RT.EQ.O.) WRITE (6,25) RC 0 131 

IF (RT.GT.G.) WRITE (6,26) RC,RT D 132 

25 FORMAT (1HL,14H FIELD WINDING// 10X19H CONDUCTOR DIAMETER, F29. 4, 1X6 D 133 

IHINCHES/) 0 134 

26 FORMAT (lHL,14h FIELD W 1 ND ING/ / lOX , 2 IH CONDUCTOR DIMENSIONS, F27. 4, D 135 

12H X,1XF6.4, 1X6HINCHES/) D 136 

WRITE (6,27) AS,PT,FE,RR,TF,FK1,PC0IL,DC0IL,BC0IL D 137 

27 FORMAT (IH ,9X15H CONDUCTOR ARE A , F 33. 4, 1X6HSQ- 1 N . / / 1 OX 13H NO. OF T D 138 

1URNS,F31.0/1CX20H MEAN LENGTH OF TURN,F27.3,2X6HINCHES//10X25H RES 0 139 
21STIVITY AT 20 DEG. C , F2 3. 4 , IX 16HM ICRO OHM INCHES/ 10X20H FIELD RES D 140 
3ISTANCE AT,F5.0,7H DEG. C, F 1 6 . 4 , 1X4H0HM S// 10X2 IH COIL INSIDE DIAME 0 141 
4TER,F26.3,2X6HINCHES/1CX22H COIL OUTSIDE D I AMETER , F25. 3/ 10X1 IH COI D 142 
5L WIDTH,F36.3) 0 143 

WRITE (6,28) DI ,DL,CL,SS,HC,SF,HV,BV,BK,WL D 144 


28 FORMAT (1H1,7H ST A TOR / / 1 CX ,2 3H STATOR INSIDE D I AMETER , F2 1 .2, 1X6HIN D 145 
1CHES/10X,24H STATOR OUTSIDE C I AMETER , F20 .2/ ICX , 32H OVERALL CORE LE D 146 
2NGTH (ONE ST ACK ) , F12 . 2/ 1 CX , 22H EFFECTIVE CORE LENGTH, F22.2/10X,17H D 147 
3 DEPTH BELOW SLOT , F2 7 . 2/ /lOX , 16H STACKING FACTOR, F28.2//10X,21H NO D 148 
4. OF COOLING DUCTS ,F 2 1 .O/lOX , 15H WIDTH OF DUCTS , F29. 2 , 1X6H INCHES// D 149 
510X.13H CORE LOSS AT,F6.1,17H KI LCL INES/SQ . IN. , F7. 1 , 2X9HWATTS/LB. ) 0 150 


IF (LTS.NE.O.) WRITE (6,29) LTS 0 151 

29 FORMAT (1CX,21H LAMINATION TH ICKNESS , F24.3 , 4H IN.) D 152 

WRITE (6,30) BP,PL,RK D 153 

30 FORMAT (1HL,6H RCTCR//10X, 16H POLE BODY WIDTH,F24.3,7H INCHES/20X, 0 154 

113H AXIAL LENGTH, F17.3/2CX,16H STACKING FACTOR ,F14. 3 ) 0 155 

IF (LTR.NE.O.) WRITE (6,31) LTR D 156 

31 FORMAT (IH ,19X,21H LAMINATION TH ICKNES S , F9 . 3 , 7H INCHES) D 157 

WRITE (6,32) PHW,PHL,RK1 0 158 

32 FORMAT (1HK,9X,16H POLE HEAD W I DTH , F24 . 3 , 7H INCHE S /20X , 1 3H AXIAL L D 159 

1ENGTH,F17. 3/20X, 16H STACKING F ACTOR , F 14. 3 ) D 160 

IF (LTRl.NE.O.) WRITE (6,31) LTRl 0 161 

WRITE (6,33) PE,HP,HP1,DR,VR,SIGMA 0 162 
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FCR^‘AT IIHK,9X,13H POLE EMBR ACE , F2 7. 3/lOX, 12H POLE HE IGHT.F28. 3, IX D 163 
16HINCHES/ lOX, 19H POLE HEIGHT ( EFF . ) , F2 I . 3/ lOX , 15H ROTOR DIAMETER, F D 164 
225.3/lOX, 17H PERIPHERAL SPEED , F20 . 0, 3X , lOH FEE T/M I N. / lOX ,23H SPEC. D 165 
3 TANGENTIAL FORCE , FI 7 . 3 , 1 1 H LBS/SQ.IN.) 0 166 

WRITE (6,34) DSH.CISH.CISHl, ALH D 167 

34 FORMAT (1HL,6H SHAFT// lOX, 28H DIAMETER (UNDER FIELD CCI 1 ) ,F13. 3, 7H D 168 

1 INCHES/10X,34H INSIDE DIAMETER (OF HOLLOW SHAFT ), F7- 3/lOX ,27H DIA D 169 
2METER (UNDER END TURN S ) , F14 . 3/ lOX, 20H LENGTH (BTW. POLES ) ,F21.3) 0 170 

IF (IBN.EG.O) WRITE (6,35) D 171 

35 FORMAT (1HL.19H CAMPER BARS (NCNE)) 0 172 

IF (DD.EQ.O..AND.IRN.NE.C) WRITE (6,36) H,B 0 173 

36 FORMAT (1HL,26H DAMPER BARS (RECTANGULAR )//lCX,22H CAMPER BAR DIME 0 174 

1NSICNS,F17.3,2H X , IX F 5 . 3 ,1X6H INCHES ) 0 175 

IF (DD.NE.O..AND.IBN.NE.C) WRITE (6,37) DD 0 176 

37 FORMAT (1HL.20H DAMPER BARS ( ROUND )// lOX ,20H DAMPER BAR DIAMETER, F 0 177 

119.3,IX6HINCHES) 0 178 

IF (IBN.NE.O) WRITE (6,38) WC ,HD , SB , TB, I BN ,RE D 179 

38 FORMAT ( IH ,9X,19H SLCT OPENING W IDTH , F20. 3/ lOX ,20H SLOT OPENING H 0 180 

1EIGHT,F19.3/10X,18H DAMPER BAR LENGTH , F21. 3/ lOX, 17H DAMPER BAR PIT D 181 
2CH,F22.3//10X,24H NO. CF DAMPER B ARS/POLE , 1 1 2// lOX ,2 5H RESISTIVITY 0 182 
3 AT 20 DEG. C,F14.3,nh MICRC-CHM INCHES) 0 183 

WRITE (6,39) TYPY D 184 

39 FORMAT (1H1,11H YOKE ( T Y PE , 1 2 , IH ) ) 0 185 

IF (TYPY-2) 40,43,48 D 186 

40 WRITE (6,41) ( ST AR ( I ) , 1= 1, 5 ) , C ASH ( 1 ) , TY , ( STAR ( 2 ) , I = 1 , 5 ) , ( DASH 1 1 ) , I D 187 

1=1,2),STAR(1) 0 188 


41 FORMAT (1HL/4(52X,1H1/),13X,5A6,5H*»» , A6 , IH-/ 13X , IH* ,31X, IH*/ 12X D 189 

1.1H*,13X,4HYCKE, 14X, 1H*,3X,4HTY =,F5.2,4H IN . / 1 IX , IH* , 30X , 2H**/l IX 0 190 
2,5A6,1H*, 1X,2A6/15X, IF*, 8X, 1H*,7X, 1H*,8X, 1H*,10X,1H1/15X,27H* STAT 0 191 
30R • FIELD • STATOR * , lOX, IH 1/ 15X , IH* , 8X ,9H* COIL * , 8X , IH* , lOX, IH D 192 
41/8X,2(7X,1H*,8X,1H* ) , lOX , IH 1 / 15X, IH* ,8X , A6 , 3H*** , 8X , IH*/ 15X , IH* , 8 0 193 


5X, 1H*,7X, 1H*,8X, 1H*/8X,2 (7X, ICH**********)///) D 194 

WRITE (6,42) CU,BCOIL 0 195 

42 FORMAT ( IHK ,9X , 20h INS I CE YOKE D I AMETER, 3X, F 7 . 3 , 7H INCHES/ lOX , 17HST 0 196 

lATCR SEPARATION, 6X,F7. 3, 7H INCHES) D 197 

GO TO 50 0 198 

43 WRITE (6,44) ( C ASH ( I ) , I = 1 , 2 ) , ( STAR ( 1 ) , I = 1 , 4 ) ,0 ASH I 1 ) 0 199 

44 FORMAT (1HL,35X,5H1 1/32X,16H 1 1 TYR/2 ( 12X , IHl ,23X, IH 1 , 0 200 

13X, IHl/) , 12X, 1H1,43X, 1H1/10X,2A6,1H-,1X,2A6,5H*****,15X, 1H1/24X, IH 0 201 
2*,15X,1H*, 15X,1H1/11X,3HTYE, 10X,1H*,5X,4HYOKE,6X,1H*,15X,1H1/2X,2( 0 202 

3I5X,A6,2H**) ,2X,A6,3H ) D 203 

45 WRITE (6,46) ( D ASH ( I ) , 1= 1 , 2 ) , ( ST AR ( I ) , 1=1 , 4 ) ,D ASH ( 1 I , ( STAR ( I ) , 1= 1 , D204 

13) D 205 

46 FORMAT (IH ,9X,2A6,5H , IX , A6 , 3H*** , 1 8X , 2HT Y/ 12X , IH 1 , 15X , IH* , 7X 0 206 

1,IH»/12X, 1H1,5X,A6,5H«**»*,1X,5HFIELD, 1X,2A6,2X,A6,3H /12X,1H1,6 0 207 

2X,lH*,8X,lH*,lX,4HCOIL,2X,lH»,aX,lH*,10X,lHl/12X,2(7X,10H* STATOR D 208 
3*) , lOX, 1H1/19X, 1H*,8X ,A6,3H***»8X, IH*, lOX, IH 1/ 12X , 2 ( 7X , 1H*8X , IH* ) , 0 209 

410X,1H1/12X,2(7X,1H*,8X,1H*)/12X,2(7X,A6,4H****)///) D 210 

WRITE (6,47) TYR,TYE, TY,CYC,CU,BCOIL 0 211 

47 FORMAT ( IHL , 9X , 3HT YR , F 30 .3 , 7H INCHES/ lOX , 3HT YE , F30. 3/lOX , 2HTY, F3 1 . D 212 

13/10X,25HDYC, YOKE INSIDE C I AME TER/ 15X, 16H ( ABO VE FLD COI L ) ,F 12. 3/1 D 213 
2CX,27HDU, STATCR OUTSIDE DI AMETER , F6. 3/ lOX , 25HBC0 I L , SPACE BTWN ST- D 214 
3AT0RS,F8.3) D 215 

GO TO 50 0 216 

48 WRITE (6,49) ( D A SH ( I ) , 1= 1 , 2 ) , ( STAR ( I ) . 1= 1, 2 ) , I DASH 1 1 ) , 1= 1 , 2 ) 0 217 

49 FORMAT ( 1 H+ , 13X , 14H ( T APEREC END S ) / IHL , 35X, IH 1 , 3X, IH 1/32X , 5H 1,3 D 218 
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1X,9HI TYR/2( 12X,ll-l,23X, IH I , 3X , IH 1/ ) , 1 2X , IH I , 43X , IH 1/ lOX , IH* ,2 

2/^6,lX,2A6,5H*»*«*, 15X,1H1/24X,1H*,15X,1H*,15X,1H1/11X,3HTYE,8X,3H# 
3**, l5X*3H***,2X,2A6t2l- — /19X , 3H*** * 8X , 4HYOKE * 9X , 3H*** / 16X , 3H**** 27 
4X,3H*** ) 

GC TC 45 

URITE (6,51) WC,RC,WI ,VsRCTCR,WYOKE,WTOTAL 

FOR^*AT (1HL»8H V^E IGHT S// lOX 1 3H STATOR COND . , El ? • 3 » 1X6HP0UNDS/ lOX 12 
IH FIELD CCND..F18.3/1CX12H STATOR I RON • F 18 • 3/ I OX , 6H R0T0R*F24,3/ 10 
2X*5H Y0KE,F25.3//10X,tF T0TAL/11X18H ( ELEC TROMAGNE T IC ) , F 1 1 . 3 ) 
fc^RITE (6,52) Cl, CP, CP, CG, 01 

FORPAT (1HL,10H CONSTANTS// lOX , 35H Cl, FUNDAMENTAL/MAX. OF FIELD F 
1LUX,F8.3/10X,18H CP, FCLE CONS T ANT , F25 . 3/ lOX , 2 ?H CM, OEM AGNET I ZA TI 
2CN FACTOR, F16.3/10X,31F CQ , CROSS MAGNETIZATION FACTOR , F 12. 3 / lOX , 2 
36H 01, POLE FACE LCSS F ACTOR , F 1 7 . 3 ) 

WRITE (6,53) AG, PC, EW, PM, P5,P6,P7 

FCRPAT (IH1,31H PERMEANCES ( L INES/AMPERE TURN ) // lOX , 8h AIR GAP,F35 
1.3,24H PER INCH OF CORE LENGTH/ lOX , 3CH WINDING LEAKAGE - STATOR SL 
2CT,F13.3/29X,10HSTAT0R END , F 14 . 3// lOX, 8H LEAKAGE/ 1 3X , 25H PM, FROM 
3R0TCR TO STAT0R/15X, 19H(BTWN. ROTOR TEETH ) , F 19 . 3/ 1 3X , 22H P5, ACROS 
4S FIELD CCIL,F18.3/13X,26H P6, FROM STATOR TC STA TOR , F 14 . 3/ 1 3X ,24H 
5 P7, STATOR TC SHAFT END,F16.3) 

WRITE (6,54) A,XR,XL,XC,XQ,XA,XB,XF,XU,SI,TC 

FCRPAT (1HL,11H R E ACT ANC ES// I0X23H AMPERE CGNDUC TORS/ INCH , F20 . 3/ 10 
1X17H REACTANCE FACTOR ♦ F26. 3// 10X23H STATOR WINDING LEAK AGE , F20 . 3 , I 
2X,7FPERCENT/10X23F ARP. REACTION (DIRECT ), F20. 3/10X22H ARM. REACTI 
3CN (QUAD. ) ,F21.3/10X21F SYNCHRONOUS ( D I RFC T ) , F22 . 3/ 10X20H SYNCHRON 
4CUS (QUAD. ) ,F23.3/10X14H FIELD LE AKAGE , F29 . 3/ 1 OX lOH TRAN S I ENT , F3 3. 
53//10X22H FIELD SELF INDLC TANCE , F2 I . 3 , IX 7HHENR 1 E S/ // 10X2 7H OPEN Cl 
6RCUIT TIME CONST ANT/ 1 7 X 1 3H (FIELD ONL Y ) , F2 3. 5 , 1X7HSECGNDS ) 

RETURN 

END 


0 219 
0 220 
0 221 
0 222 
0 223 
0 224 
0 225 
0 226 
0 227 
0 228 
0 229 
0 230 
0 231 
0 232 
0 233 
0 234 
0 235 
0 236 
0 237 
0 238 
0 239 
0 240 
0 241 
0 242 
0 243 
0 244 
0 245 
0 246 
0 247 
0 248 
0 249- 
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APPENDIX B 


DEFINITION OF FORTRAN VARIABLES 


The following is an alphabetic listing of the major FORTRAN variables used in the 
program. Tlie variables are defined and the units used in the program are given. The 
list includes approximately 75 percent of all FORTRAN variables appearing in the pro- 
gram. 


A 

AA 

AB 

ABASE 

ABODY 

AC 

ACR 

AG 

AI 

AIRGAP 

AKVA 

ALH 

ALPHAE 

ALPHAR 

ALPHAS 

ALY 

ALYC 

ALYR 

AN 

AP 

AS 

ASH 


ampere -conductors per inch of stator periphery, A/in. 
used for variety of calculations 
used for variety of calculations 

2 

area used in rotor weight calculation, in. 

2 

area used in rotor weight calculation, in. 

2 

armature conductor area, in. 

2 

effective core area per pole, in. 

specific air gap permeance per inch of core length per pole, lines/(A-turn)(in. ) 
points on material magnetization curve 
NAMELIST name 

generator output at load point G, kVA 
shaft length (between poles), in. 

temperature coefficient of resistivity of damper winding at 20° C, °C“^ 

temperature coefficient of resistivity of field winding at 20° C, °C”^ 

temperature coefficient of resistivity of armature winding at 20° C , °C 

yoke dimension used in magnetic calculation, in. 

yoke dimension used in magnetic calculation, in. 

yoke dimension used in magnetic calculation, in. 

power factor angle, rad 

2 

pole body cross-sectional area (solid area), in. 

2 

field conductor area, in. 

2 

shaft cross-sectional area, in. 
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ATH 

ATIP 

AY 

AYC 

AYR 

B 

B1 

B2 

B3 

BCL 

BCLL 

BCOIL 

BETA 

BG 

BK 

BN 

BO 

BP 

BPL 

BPLL 

BS 

BSHL 

BSHLL 

BTL 

BTLL 

BV 

BYCL 

BYCLL 

C 


2 

tooth cross-sectional area, in. 

2 

area used in rotor weight calculation, in. 
yoke area = TY*(DU+TY)*3. 14, in.^ 
yoke area = TYE*(DYC+TYE)*3. 14, in.^ 
yoke area = TYR*(DU+2. *TY)*3. 14, in. ^ 
rectangular damper bar slot width, in. 
stator slot dimension (see table Vn(c)), in. 
stator slot dimension (see table Vn(c)), in. 
stator slot dimension (see table Vn(c)), in. 

2 

core flux density at load point G, kilolines/in. 

2 

core flux density, kilollnes/in. 
field coil width, in. 

angle used in rotor weight calculations, rad 

2 

airgap flux density (no-load, rated voltage), kilolines/in. 

2 

flux density at which core loss WL is given, kilolines/in. 

number of damper bars per pole 

stator slot dimension (see table Vn(c)), in. 

pole body width, in. 

2 

pole flux density, kilolines/in. 

2 

pole flux density at load point G, kilolines/in. 

stator slot dimension (see table VII(c)), in. 

2 

shaft flux density, kilolines/in. 

2 

shaft flux density at load point G, kilolines/in. 

2 

tooth flux density at load point G, kilolines/in. 

2 

tooth flux density, kilolines/in. 
width of cooling duct, in. 

2 

yoke flux density (over field coil) at load point G, Mlolines/in. 

2 

yoke flux density, kilolines/in. 

number of parallel armature winding circuits per phase 
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Cl 

cc 

CCR 

CDD 

CE 

CF 

CK 

CL 

CM 

CONST 

CP 

CQ 

CS 

CW 

CX 

D 

D 

D 

D1 

DAMPER 

DASH 

DB 

DCOIL 

DD 

DF 

DG 

DI 

DISH 

DISHl 


ratio of fundamental' maximum to acttial maximum value of field form 
Carter coefficient (stator). 

Carter coefficient (rotor) 

2 

current density in field at load point G, A/in. 
straight portion of coil extension (see table VII(d)), in. 
pitch factor 

power factor adjustment factor 

length of one stator stack (axial direction), in. 

demagnetizing factor (direct axis) 

NAMELIST name 

ratio of average to maximum value of field form 
cross magnetizing factor (quadrature axis) 
per unit pole pitch 
winding constant 

dummy variable used in slot leakage permeance calculation 
used for distribution factor calculation 

O 

area correction for corner radii in rectangular conductor, in. 
used in interpolation between points on magnetization curve 
pole face loss factor 
NAMELIST name 

used in subroutine OUTPUT to print yoke diagram 
diameter of bender pin for forming armature coils, in. 
field coil outside diameter, in. 
damper bar diameter, in. 
distribution factor 

largest dimension of rectangular conductor (field and armature), in. 

bore (inside) diameter of stator, in. 

inside shaft diameter for hollow shaft, in. 

external shaft diameter (external to two stator stacks), in. 


DL 

DR 

DSH 

DT 

DU 

DW 

DWl 

DX 

DY 

DYC 

DZ 

E 

EB 

EC 

ED 

EDD 

EE 

EF 

EL 

EP 

ET 

EW 

EX 

EZ 

F 

FACTOR 

FCL 

FE 

FF 


damper losses at load pcmt G, W 

rotor diameter, in. 

shaft diameter (under field coil), in. 

smallest dimension of rectangular conductor (field and armature), in. 
stator outside diameter, in. 

armature winding strand diameter or width (see table Vn(d)), in. 

armature winding strand thickness (uninsulated) (see table Vn(d)) , in. 

o 

used in rectangular conductor area calculation, in. 

2 

used in rectangular conductor area calculation, in. 
yoke dimension (see table VIlQ)), in. 

2 

used in rectangular conductor area calculation, in. 
alternator efficiency at load point G, percent 
eddy factor (bottom) 

number of effective armatiire conductors 
excitation voltage at load point G, per unit 
excitation voltage, per unit 
line-to-line design voltage, rms V 
field voltage at load point G, V 
end extension length of armature coil, in. 
line -to -neutral design voltage, rms V 
eddy factor (top) 

specific stator end winding leakage permeance per inch of core length, 
lines/(A-tum)(in. ) 

eddy losses at load point G, W 

eddy factor 

frequency, Hz 

dummy variable used in slot leakage permeance calculation 

core ampere turns, A-turn 

mean length of one field coil turn, in. 

dummy variable used in slot leakage permeance calculation 



FFL 

FFLL 

FGL 

FGLL 

FGML 

FGX 

FH 

FI 

FIELD 

FKl 

FPL 

FQ 

FS 

FSC 

FSHL 

FSHLP 

FTL 

FYCL 

FYL 

FYOKE 

FYRL 

G 

GA 

GC 

GE 

GF 

GP 

GT 

GX 


total ampere turns at load point G, A-turn 

total ampere turns, A-turn 

air gap ampere turns, A-turn 

airgap ampere turns at load point G, A-turn 

demagnetization ampere turns at rated load, A-turn 

demagnetizing ampere -turns at load point G, A-turn 

airgap ampere turns (N. L. , rated volt. , for useful flux), A-turn 

field current at load point G, A 

NAMELIST name 

field winding resistance at temperature TF, ohm 
pole ampere turns, A-turn 

useful flux per pole (no-load, rated voltage), kilolines 
skew factor 

short-circuit ampere turns, A-turn 

shaft (imder field coil) ampere turns, A-turn 

shaft (xmder pole) ampere turns, A-turn 

tooth ampere turns, A-turn 

yoke ampere turns, A-turn 

yoke ampere turns, A-turn 

yoke ampere turns, A-turn 

yoke ampere turns, A-turn 

load point at which load characteristics are calculated, per unit or percent 
. 2 

airgap area, m. 

minimum air gap (air gap at center of pole) (see table Vn(e)), in. 
effective airgap, in. 

constant used in load pole-face and damper loss calculations 

maximum airgap (see table Vll(e)), in. 

ratio of slot opening width to minimum airgap 

useful flux per pole multiplying factor at load point G 
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III I Hill 


I III 



GXX 

H 

HC 

HD 

HM 

HO 

HP 

HPl 

HS 

HT 

HV 

HW 

HX 

HY 

IBN 

IDELR 

IPN 

IPX 

IQQ 

IZZ 

KSAT 

LT 

LTR 

LTRl 

LTS 

MAGNET 

OUTPUT 

P 

P5 

P6 
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flux per pole multiplying factor 

rectangular damper bar thickness, in. 

stator depth below slot, in. 

damper bar slot opening height, in. 

armature conductor length (1/2 coil length), in. 

stator slot dimension (see table VII(c)), in. 

pole height (pole body + pole head) (see table VH(g)), in. 

effective pole height, in. 

stator slot dimension (see table VH(c)), in. 

stator slot dimension (see table VII(c)), in. 

number of cooling ducts per stator stack 

stator slot dimension (see table VH(c)), in. 

stator slot dimension (see table VII(c)), in. 

stator slot dimension (see table VII(c)), in. 

number of damper bars 

voltage by which R is incremented, percent 

number of phases 

number of poles 

number of stator slots 

stator slot type 

saturation indicator (if KSAT = 0, part of alternator is saturated) 

lamination thickness (used in stacking factor calculations), in. 

pole body lamihation thickness, in. 

pole head lamination thickness, in. 

stator lamination thickness, in. 

svibroutine name 

svibroutine name 

generator input power at load point G, kW 
leakage permeance across field coil, lines/A-turn 
leakage permeance from stator to stator, lines/A-turn 



P7 leakage permeance from stator to shaft end, lines/A-turn 

PBA phase belt angle, deg 

PC specific armature slot winding leakage permeance per inch of core length, 

lines/ (A -turn) (in. ) 

PCOIL field coil inside diameter, in. 

PE pole embrace 

PF design power factor 

PH57 leakage flux across field coil, kiloline 

PH67 leakage flux from stator to stator, kiloline 

PH7L leakage flux from stator to rotor end extension, kiloline 

PHL pole head length (axial direction), in. 

PHW pole head width, in. 

PI rated line current, A 

PL pole body length (axial direction) (see table Vn(g)), in. 

PM leakage permeance from rotor to stator, lines/A-turn 

PML leakage flux from rotor to stator (see fig. 3), kiloline 

PMLA leakage flux from rotor to stator (dummy variable), kiloline 

PMLL leakage flvix at load point G, kiloline 
PP pole face losses at load point G, W 

PR field losses at load point G, W 

PS armature conductor copper losses at load point G, W 

PT nimiber of field turns 

PX number of poles 

PZ alternator losses at load point G, percent 

QAGAT air gap ampere -turns at voltage QPERV, A-turn 

QCAT core ampere-turns at voltage QPERV, A-turn 

2 

QCD flux density in core at voltage QPERV, kiloline/in. 

QFCUR field currents at voltage QPERV, A 
QN slots per pole per phase 

QPAT pole ampere -turns at voltage QPERV, A-turn 
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QPD 

2 

flux density in pole at voltage QPERV, kiloline/in. 


QPERV 

voltage at which no-load satiiration data are calculated, 

percent 

QQ 

number of slots 


QSAT 

shaft ampere-turns at voltage QPERV, A-turn 


QSD 

2 

flux density in shaft at voltage QPERV, kiloline/in. 


QTAT 

total ampere -turns at voltage QPERV, A-turn 


QTHAT 

tooth ampere -turns at voltage QPERV, A-turn 


QTHD 

2 

flux density in tooth at voltage QPERV, kiloline/in. 


QVLL 

line-to-line voltage at which no-load saturation data are calculated, rms V 

QVLN 

line -to -neutral voltage at which no-load saturation data are calculated, rms V 

QYAT 

yoke ampere -turns at voltage QPERV, A-turn 


QYD 

flux density in yoke (over field coil) at voltage QPERV, 

kiloline/in.2 

R 

alternator voltage at which no-load saturation data are calculated, per unit 

RATING 

NAMELIST name 


RC 

field coil weight, 3b 


RD 

field conductor diameter or width, in. 


RE 

damper bar resistivity at 20° C, (/Liohm)(in. ) 


RF 

type of armature winding (random or form wound) 


RGl 

armature winding resistance at temperature TST, ohm 


RK 

pole body stacldng factor 


RKl 

pole head stacldng factor 


RM 

damper bar resistivity at temperature T3 and T33, (/u. ohm) (in. ) 

ROTOR 

NAMELIST name 


RPM 

rotor rotational speed, rpm 


RR 

field coil resistivity at 20° C, (Mohm)(in. ) 


RRA 

armature winding resistance at load point G, ohm 


RRB 

field winding resistance at load point G, ohm 


RS 

armature conductor resistivity at 20° C, (Mohm)(in.) 


RT 

field conductor thickness, in. 
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J 



s 

SB 

SC 

SCR 

SD 

SF 

SH 

SHAFT 

SI 

SIGMA 

SK 

SLOTS 

SM 

SN 

SNl 

SP 

ss 

ST 

STAR 

STATE T 

STATOR 

STFK 

STRAY 

T1 

Til 

T2 

T22 

T3 

T33 


armature conductor current density at rated load, A/in.^ 

damper bar length, in. 

number of conductors per stator slot 

short circuit ratio 

distance between centerline of armature winding strands (in depth) (see 
table Vn(d)), in. 

stacking factor (stator) 

uninsulated armature winding strand height, in. 

NAMELIST name 

field self -inductance, H 

specific tangential force on rotor, psi 

stator slot skew at stator inside diameter (for one stack), in. 
NAMELIST name 

tooth width at 1/3 distance from narrowest section, in. 

strands per armature conductor in depth 

total strands per armature conductor 

total losses at load point G, W 

solid stator stack length (one stack), in. 

stator Looth losses at load point G, W 

used in subroutine OUTPUT to print yoke diagram 

number of effective armature winding turns 

NAMELIST name 

stacking factor for lamination thickness LT 
miscellaneous load losses at load point G, W 
rated-load armature temperatime, °C 
no-load armature winding temperature, °C 
rated-load field winding temperature, °C 
no-load field winding temperature, °C 
hot damper bar temperature, °C 
cold damper bar temperature, °C 
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TB 

TC 

TF 

TG 

THETA 

TP 

TS 

TST 

TT 

TTA 

TTB 

TY 

TYE 

TYPY 

SINDUC 

TYR 

UA 

VA 

VR 

WA 

WC 

WD 

WF 

WI 

WINDNG 

WL 

WN 

WO 

WPOLE 


damper bar pitch, in. 

open-circuit time constant (field only), sec 

field coil temperature at which FKl is calculated, 

total useful flux, kiloline 

angle used in rotor weight calculations, rad 

pole pitch, in. 

stator slot pitch at stator inside diameter, in. 

armature winding temperature at which RGl is calculated, 

stator slot pitch at 1/3 distance from narrow section, in. 

armature winding temperature at load point G, 

field winding temperature at load point G, *^C 

yoke dimension (see table VH(j)), in. 

yoke dimension (see table Vn(j)), in. 

type of yoke (see table VII (j)), in. 

subroutine name 

yoke dimension (see table Vn(j)), in. 

=G(M), per unit 

kilovolt-ampere rating of alternator, kVA 
rotor peripheral velocity, ft/min 
generator output power at load point G, 
stator conductor weight, lb 
no-load damper loss at temperature T3, W 
windage loss, W 
stator iron weight, lb 
NAMELIST name 

core loss at flux density BK, W/lb 
no-load pole face losses, W 
damper bar slot opening width, in. 
weight of one pole, Ib 


80 



I 


WQ 

WQL 

WROTOR 

WSHAFT 

WT 

WTOTAL 

WU 

WYOKE 

XA 

XB 

XD 

XF 

XL 

XQ 

XR 

XU 

YA 

YOKE 

YY 

ZA 

ZB 

ZC 

ZD 

ZE 

zz 

zzz 


no-load rated voltage core loss, W 
stator core losses at load point G, W 
rotor weight (=WSHAFT+PX*WPOLE), lb 
shaft weight (including portion under poles), Ib 
no-load rated voltage tooth loss, W 
total electromagnetic weight, lb 
no-load damper loss at temperature T33, W 
yoke weight, lb 

synchronous reactance (direct), percent 
synchronous reactance (quadrature), percent 
armature reaction reactance (direct), percent 
field leakage reactance, percent 
stator winding leakage reactance, percent 
armature reaction reactance (quadrature), percent 
reactance factor 

transient reactance (direct axis), percent 
=100/G 

NAMELIST name 

slots spanned per coil (number of slots between coil sides + 1) 
dummy variable used in slot leakage permeance calculation 
dummy variable used in slot leakage permeance calculation 
dummy variable used in slot leakage permeance calculation 
dummy variable used in slot leakage permeance calculation 
dummy variable used in slot leakage permeance calculation 
stator slot type (see table vn(c)) 
air gap reluctance over pole, A-turn/kiloline 
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APPENDIX C 


DEFINITION OF INPUT VARIABLES FOR EACH NAMELIST NAME 

This appendix defines all variables (FORTRAN symbols) that may be used as ii^ut to 
the homopolar inductor alternator computer program. Each variable is listed under the 
appropriate NAMELBT name. The NAMELIST names are arranged in the order in which 
the data cards must appear in the data deck. Units are given, where applicable, and each 
variable is classified as mandatory (M), conditional (C), or optional (O). A mandatory 
classification indicates that the variable must be read in. Ibe conditional classification 
indicates that, for some alternator designs, the variable is required and that, for others, 
it may be omitted. Variables identified as optional are read in at the discretion of the 
user. In each case where an optional variable is omitted, an assumption regarding that 
variable is made internal to the program. This assumption is explained in the remarks 
column of the tables. The remarks column also gives other pertinent information. 

TABLE Vn. - DEFINITIONS OF INPUT VARIABLES 
(a) NAMELIST name RATING 


mandatory; C, conditional; O, optional. 



Definition 


Kilovolt-ampere rating of alternator, kVA 

Line- to- line design voltage, rms V 

Line- to- neutral design voltage, rms V 

Frequency, Hz 

Shaft rotational speed, rpm 

Number of poles 

Design power factor 

Load points at which load characteristics 
are calculated (see sample output, 
p. 25), percent or per unit 


Classi- 

fication 

(a) 

M 

c} 

C 
C 

c 

M 

o 


Either one must be read in, or both may be 
read in 

Any two must be read in, or all three may 
be read in 


G is a subscripted variable (array size 
is 5); if not read in, program assumes 
values, 0, 0.75, 1.0, 1.25, and 1.50; 
any one or all (except 0) may be changed 
by reading in different values; program 
automatically arranges values in increas-j 
ing order; any number > 9. 0 is assumed 
to be in percent, :< 9.0 in per unit 


(b) NAMELIST name STATOR 


FORTRAN 

Definition 

Classi- 

symbol 


fication 



(a) 

DI 

Bore diameter (i. d. ), in. 

M 


DU 

Stator lamination outside diameter, in. 

M 


CL 

Length of one stator stack, in. 

M 


HV 

Number of cooling ducts 

c] 


BV 

Width of cooling duct, in. 



SF 

Stacking factor (stator) 

cl 


LTS 

Stator lamination thickness, in. 

Cj 


WL 

Core loss at flux density BK, W/lb 

M 


BK 

Flux density at which core loss WL is 

M 



given 




If there are no cooling ducts, these need 
not be read in 

Either one or both may be read in. if 
neither is read in, program assumes 
that stator in not laminated (SF = 1. 0) 


82 



TABLE Vn. - Continued. 


(c) NAMEUST name SLOTS 


HY 


1 

HX 


□ 

1 

HS 



□ 

1 



BS 




BO 

Tbs = BO) 


Types 1 and 5: Open slot, constant 
slot width. Type 5 slot is same 
as type 1, but it contains only 
one coil side. 




Type 2; Partly closed slot, 
constant slot width. 



Type 3: Partly closed slot, con- 
stant tooth width. 


Type 4: Round slot. 


FORTRAN 

[ Definition 

Classi- 

Remarks 

symbol 


! 

fication 


(a) 



(b) 


ZZ 

Slot type 


M 

See sketch 

BO 

Slot dimension, in. 1 

C 



B3 


! 





BS 







HO 







HX 







HY 







HS 


i 





HT 







IQQ 

[ Number of slots 

1 ^ J 

1 


^Variables shown in the sketch but not defined in this table are not allowable 
input. These variables are shown for reference only. 

^M, mandatory; C, conditional. 


TABLE Vn. - Continued. 


(d) NAMELIST name WINDNG 


H K 




FORTRAN 

symbol 

Definition 

Classi- 

fication 

(a) 

Remarks 

RF 

Type of coil 

M 

RF = 1. 0 for form wound coil; RF = 0 for 
random wound coil 

SC 

Number of conductors per slot 

M 


YY 

Slots spanned per coil (number of slots between 
coil sides plus one) 

M 


C 

Number of parallel circuits per phase 

M 


DW 

Strand diameter or width, in. 

M 

See sketches 

SN 

Strands per conductor in depth (radial direc- 
tion) 

C 

Read for rectangular wire only (in sketch, 
SN = 4) 

SNl 

Total strands per conductor 

M 

In sketch SNl = 8 

DWl 

Uninsulated stator strand thickness (radial di- 
rection), in. 

C 

Read for rectar^ular wire only; see 
sketches 

CE 

Straight portion of coil extension, in. 

M 

See sketches 

SD 

Distance between centerline of strands in 
depth, in. 

M 

See sketches 

PBA 

Phase belt angle, deg 

O 

j If not read in, program assumes PBA = 60^ 

SK 

Stator slot skew at stator inside diameter (for 
one stack only), in. 

O 

If not read in, program assumes SK = 0 

T1 

Rated-load armature winding temperature, ®C 

M 

Used for loss and efficiency calculations 

RS 

Armature conductor resistivity at 20° C, 
(Mohm)(in.) 

O 

If not read in, program assumes copper 
resistivity (0. 694) 

ALPHAS 

1 Armature conductor temperature coefficient 
of resistivity at 20° C , ^C 

O 

If not read in, program assumes copper 
temperature coefficient (0. 00393) 

Til 

No-load armature winding temperature, °C 

M 

Used for loss and efficiency calculations 

TST 

Armature winding temperature, °C 

O 

Program calculates and prints out arma- 
ture resistance at this temperature; if 
not read in. program assumes 
TST * 25° C 


^M, mandatory; C, conditional; O, optional. 



TABLE Vn. - Continued. 


(e) NAMELIST name ATRGAP 



(f) NAMELIST name CONST 


FORTRAN 

Definition 

Classi- 

Remarks 

symbol 


fication 




(a) 



Cl 

Ratio of fundamental maximum to actual maxi- 



Identical to those defined for convention- 


mum value of field form (field form is air- 



al salient pole alternator (ref. 7); ef- 


gap flux density distribution due to field only) 



feet of leakage flux between poles in 

CP 

Ratio of average to maximum value of field 

O 


homopolar inductor alternator is ac- 


form 


> 

counted for separately (see section 

CM 

Demagnetizing factor (direct axis) 

0 


Magnetics Calculations): if not an in- 

CQ 

Cross magnetizing factor (quadrature 

O 


put, values are calculated from for- 


axis) 



mulas given in refs. 1 and 7 

EL 

End turn length, in. 

6 ^ 


Read in if exact value is known; if not, 





program will calculate approximate 
value 

WF 

Windage losses, W ' 

O 


Read in actual value; if not read in, 





program neglects windage in effi- 
ciency calculations; to have pro- 
gram calculate approximate windage 
loss, set WF » 1.0 


^M, mandatory; C, conditional; O, optional. 
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TABLE Vn. - Continued. 


(g) NAMELIST name ROTOR 



FORTRAN 

symbol 

Definition 

Classi- 

fication 

(a) 

RK 

Stacking factor of pole body 

‘'1 

LTR 

Lamination thickness of pole body, in. 

CJ 

RKl 

Pole head stacking factor 


LTRl 

Pole head lamination thickness, in 

cJ 

PL 

Pole body length (axial direction), in. 


PHL 

Pole head length (axial direction), in. 

cJ 

PE 

Pole embrace 

""I 

PHW 

Pole head width, in. 

cJ 

BP 

Pole body width, in. 

c 

HP 

Pole height (pole body + pole head), in. 

M 

HPl 

Effective pole height, in. 

M 

WROTOR 

Rotor weight, lb 

O 

D1 

Pole face loss factor 

0 


mandatory; C, conditional; O, optional. 


Remarks 


One or both may be read in; if neither is read in program 
assumes that pole body is not laminated (RK = 1. 0) 

One or both may be read in; if neither is read in, pro- 
gram assumes solid pole head (RKl = 1. 0) 

If PL = PHL, only one (either one) need be read in; see 
sketch 

One must be read in; both may be read in 

If BP = PHW, BP need not be read in 

See sketch 

If air gap between poles is uniform, HPl = HP; if not, 

HPl > HP; Unless a better value is known, assume 
that HPl = 1. 15 HP 

If not read in, program will calculate approximate rotor 
weight 

If not read in, D1 is calculated from value of LTRl using 
the following: D1 = 1. 17 for LTRl < 0.045; D1 = 1.75 
for 0. 045 < LTRl < 0. 094; D1 = 3. 5 for 0. 094 < LTRl 
< 0. 17; D1 = 7. 0 for LTRl > 0. 17; of LTRl is not 
read in, program calculates value of LTRl based on 
RKl 



TABLE Vn. - Continued. 


(h) NAMELIST name DAMPER 


FORTRAN 

symbol 

Definition 

Classi- 

fication 

(a) 

Remarks 

BN 

Number of damper bars per pole 

M 

If BN = 0 , none of following variables for 
DAMPER need be read in 

WO 

Damper bar slot opening width, in. 

C 


HD 

Damper bar slot opening height, in. 

C 


DD 

Damper bar diameter, in. 

C 

For round damper bars only 

U 

Rectangular damper bar thickness, in. 


For rectangular damper bars only 

B 

Rectangular damper bar slot width, in. 

Cj 


SB 

Damper bar length, in. 

c 


TB 

Damper bar pitch, in. 

c 


T33 

Cold damper bar temperature, °C 

0 

If this is not read 20° C will be assumed 

T3 

Hot damper bar temperature, ®C 

c 


RE 

Damper bar resistivity at 20° C, (/j. ohm) (in. ) 

o 

0.694 Will be assumed unless otherwise 
read in 

ALPHAS 

Temperature coefficient of resistivity at 20° C, 

°c-i 

o 

0. 00393 Will be assumed unless otherwise 
read in 


^M, mandatory; C, conditional; O, optional. 


(i) NAMELIST name SHAFT 


FORTRAN 

Definition 

Classi- 

Remarks 

symbol 


fication 




(a) 


DSH 

Shaft diameter (under field coil), in. 

M 


DISH 

Inside shaft diameter (for hollow shaft), in. 

C 

Read in only for hollow shaft 

DISHl 

External shaft diameter (external to two stator 

M 



stacks), in. 



ALH 

Shaft length between poles, in. 

M 



^M, mandatory; C, conditional. 
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TABLE Vn. - Continued, 
(i) NAMELIST name YOKE 

k 


Yoke S TY 


Stator] 

Field 

coil 

Stator 


stack ' 


stack 



Type 1. 


TYR-H 


TY ? 

Field 

r—n 

T 

Stator 

coil 

stator 

1 

Inside diameter 


stack 


stack 

of yoke center 


r 


TYE 


section, DYC 


Type 2. 


* 

TYR-*- 


r 1 

r 

Yoke 

d 

TYE 

TY \ 

Field 

coil 



! 

stator 

stack 

stator 

stack 

Inside diameter 
of yoke center 


center, DYC 


Type 3: same as type 2 except that end section are 
tapered for constant flxix density. 


FORTRAN 

symbol 

Definition 

Classi- 

fication 

(a) 

Remarks 

TYPY 

Type of yoke 


M 

Three types of yokes are allowable; see sketches 

TY 

Yoke dimensions, 

in. (see sketches) 

M 


TYE 





TYR 




Needed for types 2 and 3 yokes only. 

DYC 



cj 



^M, mandatory; C, conditional. 
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TABLE VH. - Concluded. 


(k) NAMELIST name FIELD 


FORTRAN 

symbol 

Definition 

Classi- 

fication 

(a) 

Remarks 

PCOIL 

Field coil inside diameter, in. 

M 


DCOIL 

Field coil outside diameter, in. 

M 


PT 

Number of field turns 

M 


RD 

Field conductor diameter or width, in. 

M 


RT 

Field conductor thickness, in. 

C 

Do not read in for round conductors 

BCOIL 

Field coil width, in. 

C 

Do not read in if BCOIL « ALH (see 
table vn(i) 

T2 

Rated- load field ten^ierature, 

m1 

Used in loss and efficiency calculations 

T22 

No-load field temperature, °C 

Mj 

RR 

Field-coil resistivity at 20° C, (^ohm)(in. ) 

o 

If not read in, 0. 694 is assumed 

ALPHAR 

Temperature coefficient of resistivity at 20° C, 

Op-1 

o 

If not read in, 0. 00393 is assumed 

TF 

c 

Field-coil temperature, °C 

o 

Program calculates and prints out field- 
coil resistance at this temperature; if 
not read in, program assumes TF = 
25° C 


M, mandatory; C, conditional; O, optional. 
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